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ABSTRACT 

We present measurements of the stellar mass functions (SMFs) of star-forming and quiescent galaxies 
to 2; = 4 using a sample of 95 675 galaxies in the COSMOS/UltraVISTA field. Sources have been 
selected from the DRl UltraVISTA K^-band imaging which covers a unique combination of a wide 
area (1.62 deg^), to a significant depth {Ks,tot = 23.4, 90% completeness). The SMFs of the combined 
population are in good agreement with previous measurements and show that the stellar mass density 
of the universe was only 50%, 10% and 1% of its current value a.t z ^ 1.0, 2.0, and 3.5, respectively. 
The quiescent population drives most of the overall growth, with the stellar mass density of these 
galaxies increasing by 0.64^0 03, 1.461q q7, 2.7llg 22 dex since z = 1.0, 2.0, and 3.5, respectively. 

By contrast, the mass density of star-forming galaxies grows by 0.40lg'Q3, 0.87^^06' ^^'^ 
dex over the same redshift range. At z > 2.5, star-forming galaxies dominate the total SMF at all 
stellar masses, although a nonzero population of quiescent galaxies persists to z = 4. Comparisons 
of the Ks-selected star- forming galaxy SMFs to UV-selected SMFs at 2.5 < z < 4 show reasonable 
agreement and suggests UV-selected samples are representative of the majority of the stellar mass 
density at z > 3.5. We estimate the average mass growth of individual galaxies by selecting galaxies 
at fixed cumulative number density. The average galaxy with Log(M*/MQ) = 11.5 at 2: = 0.3 has 
grown in mass by only 0.2 dex (0.3 dex) since z — 2.0(3.5), whereas those with Log(M*/M0) — 
10.5 have grown by > 1.0 dex since z = 2. At z < 2, the time derivatives of the mass growth are 
always larger for lower-mass galaxies, which demonstrates that the mass growth in galaxies since 
that redshift is mass-dependent and primarily bottom-up. Lastly, we examine potential sources of 
systematic uncertainties on the SMFs and find that those from photo-z templates, SPS modeling, and 
the definition of quiescent galaxies dominate the total error budget in the SMFs. 
Subject headings: galaxies: mass function - galaxies: evolution - galaxies: high-redshift - galaxies: 
fundamental parameters 



1. INTRODUCTION 

In the current ACDM paradigm, the dominant 
structures in the universe are dark matter halos which 
grow out of an initial field of density perturbations via 
gravitational collapse (Whi te fc Rees,1978. ) . Simulations 
and analytical models show that this process proceeds 
primarily in a hierarchical, bottom-up manner, with 
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low-mass halos forming early and subsequently growing 
via continued accretion and m erging to form more 
massive halos at later t imes (iWhite fc FrenkI 119911 : 
IKauffmann fc Whitd [19931 : IKauffmann et al.lll999n . 

In contrast to the predicted hierarchical growth of the 
dark matter halos, observational studies suggest that the 
stellar baryonic component of the halos (i.e., galaxies), 
may grow in an anti-hierarchical, top-down manner. 
It appears that many of the most massive galaxies 
(Log{M star /^q) > 11) in the local universe assembled 
their stellar mass rapidly and at early times (z > 
2), whereas lower-mass galaxies grew m ore gradually 
over cosmic time (e.g.. iMarchesini et al . 2009, 2010t 
Ilbert et all 120101: ICaputi et all 120111 : IBranimer et ahl 
20TW . 

Understanding these apparently contrasting evolution- 
ary paths for the dark matter assembly and stellar mass 
assembly of galaxies is a significant challe nge for current 
models of galaxy form ation (e.g., Ma rchesini et al.ll2009l : 
iFontanot et al.l [20091 ) . In particular, the differential 
evolution between the baryonic and non-baryonic com- 
ponents of galaxies makes it clear that the baryonic 
physics of galaxy formation must be more complex 
than the cooling of gas onto halos at a rate dictated by 
gravity. Indeed, it implies that there is a tenuous balance 
between gas accretion rates, gas consumption rates (in 
both star formation events and black hole growth). 
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mergers, as well as feedback processes such as AGN 
activity, supernovae, or stellar winds. It is also clear 
that the efficiency of these processes mu st scale with 
halo mass and evolve with re d shift Ce.g.. iSchaye et al.l 
[Mot iWeinmann elall [20T1 iHenriques et al.l I2012D . 
Given the complex, non-linear interplay between these 
processes, and the various possible prescriptions of 
implementing them within models, it is important to 
have a benchmark for the models so that we can evaluate 
if progress is being made. 

For cosmological simulations, the benchmark that has 
been most widely adopted is the ability of models to 
match the volume density of galaxies as function of their 
stellar mass, also known as the stellar mass function 
(hereafter SMF). If a model can reproduce the SMFs at 
various redshifts it suggests (although does not prove) 
that it may be a better description of the baryonic 
physics of galaxy formation than those that do not. 
Given that it is a key benchmark for models, the most 
precise and accurate measurements of the SMF possible 
over as large a range in redshift and stellar mass are 
valuable quantities. 

In recent years, with the growth of deep and 
wide-field near infrared (NIR) imaging surveys, there 
have been myriad measurements of the e volution of 
the S MFs from the local univer s e (e. g .. iCole et al 
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Bundv et a l. 2006; Pozzctti ct al. 2007; Arnouts ct al 
2007t iPerez-Gonzalez et al .2008.: Drorv et al. 2009 : 
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Bielbv et alj |2012l: iMoustakas ct al. 20131: lllbert et alJ 
20131 ). In this paper we present an improved mea- 
surement of the SMF of galaxies over the redshift 
range 0.2 < z < 4.0. These measurements are made 
from a new Kg-selected catalog of the COSMOS field 
whic h uses data from the DRl UltraVISTA survey 
fsee iMcCracken et al.|[20Tl ). The UltraVISTA catalog 
is unique in its combination of covering a wide area 
(1.62 deg^), to a relatively deep depth (Ks.tot < 23.4, 
90% completeness). This combination allows the most 
accurate measurement of the high mass end of the SMFs 
up to z = 4.0 to date. Details of the catalog and a 
public release of all catalog data products are presented 
in a companion paper by Muzzin et al. (2013). 

We note that an independent analysis of the SMFs out 
to z = 4 using the UltraVISTA data has also recently 
been performed by lllbert et"al] ()2013[ ). That analysis 
is based on a different catalog than that of the Muzzin 
et al. (2013), and uses different photometric redshift 
and stellar mass fitting techniques. In an appendix we 
make a more detail ed comparison bet ween our SMFs 
and those derived by lllbert et~al\ (|2013t ) . 

The layout of this paper is as follows. In § 2 we 
present details of the COSMOS/UltraVISTA dataset 
and discuss the stellar mass and photometric redshift 
measurements. In § 3 we detail how the SMFs and the 
uncertainties are calculated. In § 4 we derive the SMFs 
of star forming and quiescent galaxies and the stellar 
mass density and number density evolution up to z = 4. 
In § 5 we present a discussion of our results, including 
a comparison to UV-selected SMFs at z > 3 and an 
estimation of the typical mass growth of galaxies using a 



fixed cumulative number density approach. We conclude 
in § 6 with a summary of our results. In an appendix we 
present a detailed look at possible sources of systematic 
error and their effect on the derived SMFs. Throughout 
this paper we assume a JIa = 0.7, fi™ = 0.3, and Hq = 
70 km s~^ Mpc~^ cosmology. All magnitudes are in the 
AB system. 

2. THE DATASET 

This study is based on a Ks-selected catalog of 
the COSMOS/UhraVISTA field from Muzzin et al. 
(2013). The catalog contains PSF-matched photometry 
in 30 photometric bands covering the wavelength range 
O.lSAtm — >■ 24Mm a nd includes the available GALEX 
(IMartin et al.l[2005l). CFH T /Su baru llCapak et al.ll2007D , 



UltraVISTA (IMcCracken et al.l |2012D, and S^SMOS 
(^San ders et al.l I2007D datasets. Sources are se- 
lected from the DRl UltraVISTA K^-band imaging 
(jMcCracken et al.[l2012t ) which reaches a depth of Ks^tot 
< 23.4 at 90% completeness. A detailed description of 
the photometric catalog construction, photometric red- 
shift (zphot) measurements, and stellar mass (hereafter, 
Mstar) estimates is presented in Muzzin et al. (2013). 
A public release of all data products from the catalog is 
also presented with that paper. Here we briefiy describe 
the aspects of the catalog relevant to the measurement 
of the SMFs. 

2.1. Photometric Redshifts and Stellar Masses 

Each galaxy in the catalog has a Zphot determined by 
fitting the photometry in the 0.15/im — !■ 8.0/im bands 
to template Spectral Energy Distributions ( SEDs) us- 
ing the EAZY code ([Brammer et al.l 120081) . In de- 
fault mode, EAZY fits photometric redshifts using linear 
combinations of 6 templa tes from the PEGASE models 
(jFioc fc Rocca-Volmeran gc 1999) as well as an additional 
red template from the [Marastorj ()2005D models. In or- 
der to improve the accuracy of the Zphot for high-redshift 
galaxies we added two new templates to the default set, a 
~ 1 Gyr old poststarburst template, as well as a slightly 
dust-reddened Lyman Break template (see Muzzin et al. 
2013). Comparison of the Zphot to 5100 spectr oscopic red- 
shifts from the zCOSMOS-bright 10k sample ([Lilly et al.l 
|2007[) , as well as 19 spectros copic r edshifts for red galax- 
ies at z > 1 (Ivan d c Sand e et al.l [201 1: Onodcr a et al.l 
I2012t IBezanson et"al.l 120131: Ivan deSande et al.l l2012| l 
shows that the Zphot have an rms dispersion of ^z/(l 
+ z) — 0.013 and a > 3(t catastrophic outlier fraction of 
1.6%. 

Stellar masses for all galaxies have been determined 
by fitting the SEDs of galaxies to stellar popula- 
tion synth esis (SPS) models using the FAST code 
([Kriek et a l. 2009). It is well-known that the M^tar de- 
rived from SED fitting depends on the assumptions made 
(met allici ty, SPS model, dust law, IMF) in this pr o- 
cess (e.g.. [Marchesini et al.ll2009l: iMuzzin et al.ll2009al[bl : 
iConrov et al.ll2009D . These assumptions typically result 
in systematic chang es to the SMFs, rather than larger 
random errors (e.g.. [Marchesini et al.l 120091 ). Given the 
complexity of these systematic dependencies, in this pa- 
per we base the majority of the analysis on a default set 
of assumptions for the SED modeling and then in the 
Appendix we expand the range of SED modeling param- 
eter space and explore their effects on the SMFs. 
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Fig. 1. — Example SEDs of red and blue galaxies in three redshift ranges: 2.5 < z < 3.0 (top row), 3.0 < z < 3.5 (middle row), 3.5 

< 2 < 4.0 (bottom row). The second and fourth columns show galaxies that have magnitudes near the limiting magnitude of the SMFs 
(Ks ~ 23.4), and the first and third columns show galaxies that are ~ 1 mag brighter. There are no red galaxies with Ks ~ 22.4 and 3.5 

< 2 < 4.0. The SEDs of galaxies at Ks = 23.4 limit typical have S/N ~ 5 and therefore we have limited the SMFs to a limiting Mstar that 
corresponds to this limit. 



For the default set of Mstar we fit the SEDs to a 
set of models with exponentially-declining star forma- 
tion histories of the form SFR cx e"*/"^, where t is the 
time since the onset of star formation, and r sets the 
ti mescale of the decline in t he SFR. We use the models 
of lBruzualfc C hariot' (2003), hereafter BC03, with solar 
metallicitv, alCalzctti ct al. (2000) dust law, and assume 
a IKroupal (|2001 ;) IMlQ. We allow log(T/Gyr) to range 
between 7.0 and 10.0 Gyr, log(i/Gyr) between 7.0 and 
10.1 Gyr, and A^, between and 4. The maximum al- 
lowed age of galaxies is set by the age of the universe at 
their Zphot ■ Further details on the default model set and 
the fitting process are discussed in Muzzin et al. (2013). 

3. CONSTRUCTION OF THE STELLAR MASS FUNCTIONS 

Here we outline how the SMFs for the quiescent, star 
forming, and combined populations are constructed. 

3.1. Galaxy Sample and Completeness 

The Muzzin et al. (2013) catalog contains a total of 
262 615 objects down to a 3(t limit of Kg < 24.35 in 
a 2.1" aperture. From that parent sample we define a 

The stellar m asses in the Muzzin et al. (2013) catalog are 
computed with a Chabrier (2003) IMF. Fo r easy compari son with 
the literature we have converted them to a IKroupal l|20011 ) IMF by 
increasing them by 0.04 dex 
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Fig. 2. — Grayscale representation of the density of galaxy stel- 
lar masses as a function of redshift in the Ks-selected catalog. 
The 100% and 95% mass-completeness limits determined using 
the deeper datasets are shown as the purple and red curves, re- 
spectively. Also shown is the 100% completeness limit for an SSP 
formed at 2 = 10. 

mass-complete sample for computing the SMFs by ap- 
plying various cuts to the catalog. 

Simulations of the catalog completeness (see Muzzin 
et al. 2013, Figure 4), show that the 90% point-source 
completeness limit in total magnitudes for the UltraV- 
ISTA data is Kg.tot = 23.4 after the blending of sources 
is accounted for. This limit in Ks^tot also corresponds 
to the ~ 5a limit for the photometry in the 2.1" color 
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Fig. 3. — UVJ color-color diagram at various redshifts for galaxies more massive than the 95% mass-completeness limits. The bimodality 
in the galaxy population is clearly visible \ip to z = 2. The cuts used to separate star forming from quiescent galaxies for the SMFs are 
shown as the solid lines. 



aperture, and therefore is a sensible limiting magnitude 
for computing the SMFs. 

As a demonstration of the quality of the SEDs near 
the 90% completeness limit, in Figure 1 we plot some 
randomly-chosen examples of red and blue galaxy SEDs 
in three redshift bins: 2.5 < z < 3.0 (top row), 3.0 < 2 < 
3.5 (middle row), 3.5 < z < 4.0 (bottom row). We plot 
SEDs of galaxies that have fluxes near the 90% complete- 
ness limit {Ks^tot ~ 23.4), as well as SEDs of galaxies 
that are ~ 1 magnitude brighter {Ks.tot ~ 22.4). Figure 
1 shows that the SEDs of both red and blue galaxies at 
Ks,toi ^ 22.4 are very well constrained. It also shows 
that at Ks^tot ~ 23.4, the SEDs are also reasonably well- 
constrained; however, the typical S/N in a 2.1" aperture 
is ~ 5. 

It is possible to include galaxies fainter than the 90% 
i^s,tot completeness limit in the SMFs and correct for this 
incompleteness; however, given that the quality of the 
SEDs near Kg^tot ~ 23.4 becomes marginal, we have cho- 
sen to restrict the sample to galaxies with good S /N pho- 
tometry. This ensures that all galaxies included in the 
SMFs have reasonable well-determined M^tar and Zphot- 

When constructing the SMFs we also exclude objects 
flagged as stars (star = 1) based on a color-color cut, as 
well as those with badly contaminated photometry (SEx- 
tractor flag K_flag > 4). Objects nearby very bright 
stars (contamination — 1) or bad regions (nan_contain 
> 3) are also excluded, and the reduction in area from 
these effects is taken account of in the total survey vol- 
ume. 

Once these cuts are applied, the final sample of galaxies 



available for the analysis is 160 070. In Figure 2 we plot 
a grayscale representation of the yistar of this sample as 
a function of Zphot- In general, the sample is dominated 
by objects at z < 2; however, there are reliable sources 
out to z = 4. 



3.2. Stellar Mass Completeness vs. z 

Figure 2 shows the Mgtar of galaxies down to 90% Kg- 
band completeness limit of the survey; however, in order 
to construct the SMFs, the limiting Mstar above which 
the magnitude-limited sample is complete needs to be 
determined. In order to estimate the redshift-dependent 
comple teness limit in M,;j„^ we adopt the approach devel- 
oped in lMarchesini et al.l (|2009( ). which exploits the avail- 
ability of other survey data that are deeper than Ultra- 
yiSTA. Specifically, we employed the JC-selected FIRES 
(|Labbc ct al."2003':'Forster Schre iber et al.ll2006D and the 
FIREWORKS ( Wuyts et al. 200$ catalogs, aheady used 
in Marchc sini et al.l (|200^ 12010 !). and the iJigo-selected 
catalo gs over the Hubble U ltra-Deep Field (H-UDF) used 
in M archesini et all (|20ll) . The FIRES-HDFS, FIRES- 
MS1054, FIREWORKS, and HUDF reach limiting mag- 
nitudes of i^s^tot = 25.6, 24.1, 23.7, and 25.6, respec- 
tively. 

Briefly, to estimate the redshift-dependent stellar mass 
completeness limit of the UltraVISTA sample at Ks,tot = 
23.4, we first selected galaxies belonging to the available 
deeper samples. We then scaled their fiuxes and Mgtar to 
match the i^-band completeness limit of the UltraVISTA 
sample. The upper envelope of points in the {Mstar,scaied 
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- z) space, encompassing 100% of the points, represents 
the most massive galaxies at = 23.4, and so pro- 
vides a redshift-dependent Mgtar com pleteness limit for 
the U ltraVISTA sample. We refer to iMarchesini et al.l 
^00§) for a more detailed description of this method. 
In Figure 2 we show this empirically-derived 100% mass- 
completeness limit as the purple curve. Also, for refer- 
ence we show the mass-completeness limit for a simple 
stellar population (SSP) formed at z = 10 which is ex- 
treme but indicative of a maximally-old population. 

Figure 2 shows that for galaxies at z < 1.5 and Kg^tot 
~ 23.4, the most extreme M/L ratios are less extreme 
than an SSP. Around z = 1.5 the SSP curve and the em- 
pirical 100% completeness curve cross each other which 
implies that there exist galaxies that have larger M/L 
ratios than an SSP. Such galaxies are typically galaxies 
with intermediate-to-old ages (for their redshift) with up 
to several magnitudes of dust extinction. More detailed 
SED modeling for galaxies with spectroscopic redshifts 
shows that these dusty and old galaxies are not uncom- 
mon among the mas sive galaxy population at 2; > 1.5 
(see e.g., IKriek et al.ll2006[ I2OO8I : iMuzzin et ani2009al ). 

As Figure 2 shows, the empirically-derived 100% mass- 
completeness limits are high due to the old and dusty 
population. Adopting the empirical 100% completeness 
limit for the SMFs therefore requires the exclusion of 
57% of the magnitude-limited sample. 

The 100% completeness limit is set by most extreme 
M/L ratio at any given redshift, regardless of the fre- 
quency of its occurrence. In principle, if only a small 
fraction of objects have these extreme M/L ratios, then 
adopting the 100% mass-completeness limit is an ineffi- 
cient use of the data. In order to try to make better use 
of the dataset we also derived 95% mass-completeness 
limits for the sample and this limit is also plotted in Fig- 
ure 2. 

At all redshifts the 95% mass-completeness limits are 
0.2-0.3 dex lower showing that it is only a small fraction 
of the overall population of galaxies that have extreme 
M/L ratios. If we adopt the 100% mass-completeness 
limits, the resulting sample of galaxies is 67 942. Adopt- 
ing the 95% mass-completeness limits increases the sam- 
ple by a factor of 1.4 to 95 675 galaxies. Given this sub- 
stantial increase in statistics, and the advantage gained 
by probing further down the SMFs at higher redshift, we 
have adopted the 95% mass-completeness limits for the 
SMFs, but correct the lowest-mass bin in each SMF by 
5% in order to account for this. 

3.3. Separation of Quiescent and Star-Forming Galaxies 

It is well-known that the overall galaxy popula- 
tion is bi- modal in the distribut i on of color s and 
SFRs (e.g ., Kauffmann et al. 2003: 'Hogg et aL| [20p: 
iBalogh et al.,.2004: .Blanton &: Moustakas,2009i rand that 
this bi-modality persists out to high redshift (e.g. . 
Bell et al."2004'; 'Taylor gtal.][2009l: iWilliams et al|[2009l: 
Brammcr et al. 2009, 20 lit). Given the bi-modality, sep- 
arating the evolution of the SMFs of star-forming and 
quiescent galaxies as a function of redshift is useful for 
understanding the relationship between the two popula- 
tions. 

In recent years, several methods have been developed 
to classify galaxies into these categories. In this anal- 
ysis we perform classification between the types using 



the rest-frame U - V vs. V - J color-color diagram (here- 
after the UV J diagram) . The UV J cla ssification has been 
used in many previous studies ( e.g., 'Labbe e t al.l 120051: 
Wuvts et al.l 120071 : IWilliams et al . 2009; Bra mmer et~al] 
20111: iPatel et al.ll20ll ^ The main advantage of this 
method, as opposed to a cut in a derived quantity such 
as specific star formation rate (SSFR), is that it is clos- 
est to a direct observable, namely the color of a galaxy. 
SSFRs derived from SED fitting depend strongly on the 
chosen form of the star formation history (for both the 
Mstar and SFR) (e.g., [Maraston et al. 2010) and SFRs 
derived from Ha, or mid-IR flux can be contaminated 
by AGN activity. Rest-frame colors do not require SED 
fitting, and AGN tend not to contaminate these bands 
significant unless they dominate the bolometric luminos- 
ity. 

In Figure 3, we plot the U - V vs. V - J diagram for 
galaxies more massive than the 95% mass-completeness 
limits in several redshift bins. The galaxy bi-modality 
is clearly visible in the UVJ diagram up to z = 2, but 
thereafter becomes less pronounced at the Mstar com- 
pleteness limits probed by the Kj-selected UltraVISTA 
catalog. 

To distinguish between star forming and quiescent 
galaxies we use box regions in the UVJ diagram that 
are similar, altho u gh n ot identical to those d efined 
in 'Williams et a.1.1 (I2009D : iWhitaker et~al] ()201lD and 
[Brammcr et al. (j201ll ). These regions are plotted as the 
solid lines in Figure 3. Quiescent galaxies are defined as, 

U -V > 1.3, V - J < 1.5, [all redshifts] (1) 
U -V > {V -J)x 0.88 + 0.69, [0.0 <z < 1.0] (2) 
U-V > [V ~J)x 0.88 -I- 0.59, [1.0 <z < 4.0] (3) 

We note that these boxes are chosen arbitrarily, with 
the main criteria being that they lie roughly between the 
two modes of the pop ulation seen in Figure 3. They were 
originally defined by IWilliams et al.l ()2009D who defined 
them in such a way to maximize the difference in SSFRs 
between the regions; however, ou r rest-frame color dis - 
tribution is slightly different than IWilliams et al.l ()2009f ) 
which is the reason that we have adjusted the box lo- 
cations. In the appendix we explore the effect on the 
SMFs of moving the location of the boxes in UVJ space. 
In general, we find that changing the UVJ box has little 
effect on the high-mass end of the quiescent SMF, and 
the low-mass end of the star-forming SMF because those 
galaxies are dominated by very red and very blue galax- 
ies, respectively. It docs have a larger effect on the SMF 
for intermediate mass galaxies, which is not unexpected 
given that these are typically the transition population 
at most redshifts. 

3.4. Stellar Mass Function Construction and Fitting 

With well-defined mass-completeness limits as a func- 
tion of redshift and criteria for separating star-forming 
and quiescent galaxies, SMFs can now be computed. 
We employ two methods to determine the SMFs, the 
l/Vmax method and a maximum-likelihood method. 
These methods have different strengths and weaknesses. 
The 1/Vmax method has the advantage that it does not 
assume a parametric form of the SMF, allowing a direct 
visualization of the data; however, it is a fully normal- 
ized solution and is susceptible to the effects of clustering. 
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Fig. 4. — Uncertainty in the number density of galaxies with 
Log(Mstar/M0) = 11.0 due to cosmic varia nce as as a fun c tion o f 
redshift calculated using the prescription of IMoster et al.l 11201 II ). 
Other surveys with smaller areas but also more independent sight 
lines are shown for comparison (see text for details). The uncer- 
tainties in UltraVISTA due to cosmic variance are ~ 8 - 15% at 
Log(Mstar/M0) = 11.0 over the full redshift range. 

Conversely, the maximum-likelihood method has the ad- 
vantage that it is not affected by density inhomogeneities 
(e.g., Efstathiou ct al. 1988); however, it does assume a 
functional form for the fit. 

3.4.1. The 1/Vmax Method 

To measure the SMFs for the sample we have applied 
an ex tended versio n of the l/Yraax algor ithm (Schmidl 
[19681) as defined in lAvni fc Bahcaiil ([l980[ ). The method 
has been used to determine the rest- frame optica l lumi- 
nosity func t ions and SMFs by iMarchesini et al.l (|2007L 
[2(M [201(1 [20ll) and we refer to those papers for an 
in-detail description of the method. 

In brief, for each yistar we determine the maximum 
volume within which an object of that yistar could be de- 
tected. This volume is determined as a function of M.star 
using the maximum redshift that the survey is complete 
for objects of that M^tor- The SMF is then calculated by 
counting galaxies in bins of yistar and correcting those 
bins with 1/Vmax- Poisson error bars are determ ined for 
each bin using the prescription of iGehreTsI (|1986[ ). which 
is valid for small number statistics. 

3.4.2. The Maximum Likelihood Method 

The SMFs are also determ ined using the maximum- 
likelihood method outlined bv lSandage et al.l ()1979D . For 
this method it is assumed that the n umber density o f 
galaxies {^{Mstar)) is described by a iSchechteil (|1976l ) 
function of the form, 

$(A/) = (/nlO)$*[10(*^-*'^*)(i+")] X ea;p[-10(*^-^^*)], 

(4) 

where M = \og{Mstar/^Q), a is the low-mass-end slope, 
M* = log(M*j^,,/Mo) is the characteristic mass, and $* 
is the normalization. For each possible combination of a 
and y^star likelihood that each galaxy would be found 
in the survey is calculated. The best-fit solution for a 



and y^*star i'^ each redshift bin is obtained by maximizing 
the combined likelihoods of all galaxies (A) with respect 
to these parameters. The is determined by requiring 
that the total number of observed galaxies is reproduced. 
The errors in $* are then determined from the minimum 
and maximum values of $* allowed by the confidence 
contours in the a vs. M*^^^ pla ne. Further details of the 
fittin g process can be found in IMarchesini et al.l ()2007l 
[20091) . 

3.4.3. The Low-Mass-End Slope a 

The SMFs are computed over a large redshift range, 
and as was shown in Figure 2, the limiting 95% com- 
pleteness limit in ^Istar is a strong function of redshift. 
The SMFs reach ~ 1.5 dex deeper than M*i^^ at z < 0.5, 
but only to ^ M*^^^ itself at z = 3.5. This means that 
a is well constrained at z < 2, but is poorly constrained 
at z > 2. Given the well-known correlation between a 
and M*(jj^, it is important to be aware that the true 
uncertainties in quantities such as M*^^^, or the stellar 
mass density can be systematically larger than the ran- 
dom uncertainties due to the data not reaching a mass 
limit sufficiently low to constrain a. In order to quan- 
tify the uncertainties in all parameters better, in all red- 
shift ranges we have performed the Schechter function 
fits both with a as a free parameter, and fixing it to a 
known value. 

In the fits with a held fixed, we have chosen values of 
a = —1.2, —0.4, and —1.3 for the total, quiescent, and 
star forming populations, respectively. As discussed in § 
4, these values are similar to those derived at z < 1 when 
a is fit as a free parameter. They are also consistent with 
values at z > 1 derived from studies that probe the low- 
mass end better than UltraVISTA (e.g., Fontana et aU 
2006; Perez -Gonza lez et al . 2008; Marcliesini et al. 200^ 
Stark et al.ll2009HLee et alJ l2012'). 

In addition to fits with a fixed and free, we have also 
performed fits to a "double" Schechter function. Several 
recent studies have shown that at low- redshift, the low- 
mass end of the SMF {ljOg{Mstarly^Q) < 9.5) is better 
described by the sum of two Schechter fu nctions with 
ident i cal M*,„^, but diffe rent $* and a (e.g.. lLi fc White! 
[2009l:IBaldrv et al.l[2QT2[) . The UhraVISTA SMFs reach 
the limiting M^^tar where a clear departure from a single 
Schechter function fit is seen at z < 1, and so for the 
SMFs in that redshift range fits are also performed with 
a double Schechter function. 

3.4.4. Determination of Uncertainties in the SMFs 

In addition to the Poisson uncertainties, there are sev- 
eral other sources of uncertainty in the construction of 
SMFs that need to be taken into account. The largest 
are caused by the fact that M-star itself is not an ob- 
servable quantity, but is derived from observables (i.e., 
multiwavelength photometry) using a set of models. The 
effect of photometric uncertainties on the derived Zphot 
and yistar is a non-trivial function of color, magnitude, 
and redshift caused by a range of data depths in various 
bands within the survey. These types of uncertainties 
can be substantial for different galaxy types and have 
been typically neglected in most previous measurements 
of the SMFs. 

In order to calculate uncertainties in the SMFs due to 
photometric uncertainties we perform 100 Monte Carlo 
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Fig. 5. — Stellar mass functions of all galaxies, quiescent galaxies, and star-forming galaxies in different redshift intervals. The 
shaded/hatched regions represent the total Icr uncertainties of the maximum-likelihood analysis, including cosmic variance and the er- 
rors from photometric uncertainties as derived using the MC realizations. The normalization of the SMF of quiescent galaxies evolves 
rapidly with redshift, whereas the normalization for star-forming galaxies evolves relatively slowly. In particular, there is almost no change 
at the high-mass end of the star forming SMF, whereas there is clear growth at the high-mass end of the quiescent population. There is 
also evidence for evolution of the low-mass end slope for quiescent galaxies. At low-redshift a double Schechter function fit is required to 
reproduce the total SMF. 

metric uncertainties as derived using the MC realiza- 
tions, and the error due to cosmic variance. For the 
maximum-hkehhood method, the total la random errors 
of the Schechter function parameters a, M*^^^,, and $* 
are the quadrature sum of the errors from the maximum- 
likelihood analysis, the errors from photometric uncer- 
tainties as derived using the MC realizations, and the 
error due to cosmic variance (affecting only the normal- 
ization $*). 



(MC) realizations of the catalog. Within each realiza- 
tion the photometry in the catalog is perturbed using 
the measured photometric uncertainties. New Zphot and 
Mstar are calculated for each galaxy using the perturbed 
catalog. The 100 MC catalogs are then used to recalcu- 
late the SMFs and the range of values gives an empirical 
estimate of the uncertainties in the SMFs due to un- 
certainties in Mstar and Zphot that propagate from the 
photometric uncertainties. 

In addition to these Zphot and M^tar uncertainties, the 
uncertainty from cosmi c variance is also included us- 
ing the prescriptions of iMoster et al.l ()2011l ). In Figure 
4 we plot the uncertainty in the abundance of galax- 
ies with Log(Mstcir/M0) = 11.0 due to cosmic variance 
as a function of redshift. Cosmic variance is most pro- 
nounced at the high-mass end where galaxies are more 
clustered, and at low redshift, where the survey volume 
is smallest. Also plotted in Figure 4 are the cosmic vari- 
ance uncertainties from other NIR surveys such as FIRE 
WORKS 
200 



(Wuvtsctal. 20Qa), MUYSC (Quad ri et al 
iMarchesini et all 1200911 ■ NMBS (fWhitaker et al 



20m. and the UPS ( iWiUiams et al.l l2009D. These sur 
veys cover areas that are factor of ~ 50, 16, 4, and 2 
smaller than UltraVISTA, respectively. Figure 4 shows 
that the improved area from UltraVISTA offers a factor 
of 1.5 improvement in the uncertainties in cosmic vari- 
ance compared to even the best previous surveys, and 
that over the full redshift range the uncertainty from 
cosmic variance is ~ 8 - 15% at Log(Mstar/M0) — 11.0. 

The total uncertainties in the determination of the 
SMFs are derived as follows. For the 1/Vmax method, 
the total la random error in each mass bin is the quadra- 
ture sum of the Poisson error, the error from photo- 



4. THE STELLAR MASS FUNCTIONS, MASS DENSITIES 
AND NUMBER DENSITIES TO Z = 4 

4.1. The Stellar Mass Functions 

In Figure 5 we plot the best-fit maximum-likelihood 
SMFs for the star-forming, quiescent, and combined pop- 
ulations of galaxies. Figure 5 illustrates the redshift evo- 
lution of the SMFs of the individual populations, which 
we discuss in detail in § 5. To better illustrate the relative 
contribution of both star-forming and quiescent galaxies 
to the combined SMF, in Figure 6 we plot the SMFs de- 
rived using the l/Vmax method (points), as well as the 
fits from the maximum-likelihood method (filled regions) 
in the same redshift bins. The SMFs of the combined 
population are plotted in the top panels, and the SMFs 
of the star- forming and quiescent populations are plotted 
in the middle panels. Within each of the higher redshift 
bins, the SMFs from the lowest-redshift bin (0.2 < z < 
0.5) are shown as the dotted line as a fiducial to demon- 
strate the relative evolution of the SMFs. The fraction 
of quiescent galaxies as a function of Mstar is shown in 
the bottom panels and the best-fit Schechter function pa- 
rameters for these redshift ranges are listed in Table 1. 

For reference, in the lowest-redshift panel (0.2 < z < 
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Fig. 6. — Top panels: The SMFs of all galaxies in different redsliift bins from 0.2 < z < 4.0. The black points represent the SMFs 
determined using the l/Ymax method and the black solid curves are the SMFs determined using the maximum-likelihood method. The 
gray shaded regions represent the total Icr uncertainties of the maximum-likelihood analysis, including cosmic variance and the errors from 
photometric uncertainties as derived from the MC simulations. Overplotted in the 0.2 < z < 0.5 bin are the SMFs from Cole et al. (2001), 
Bell et al. (2003), and Baldry et al. (2012). In the remaining redshift bins the dotted curve is the total SMF from Ultra VISTA in the 0.2 
< z < 0.5 bin. Middle panels: SMFs as with the top panels, but for the quiescent galaxies (red points, red solid curves) and star-forming 
galaxies (blue points, blue solid curve). The orange and cyan shaded regions represent the total Icr uncertainties of the maximum-likelihood 
analysis for quiescent and star-forming galaxies, respectively. Bottom panels: Fraction of quiescent galaxies as a function of M.gtar- 
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Fig. 7. — Left panel: The evolution of the stellar mass density of galaxies from z = 4 to z = down to a limit of log(Mstar /M.q) = 8.0. 
The UltraVISTA measurements are shown in red with error bars representing total la random errors inclusive of cosmic variance and the 
errors from photometric uncertainties as derived using the MC simulations. Other measurements from the literature are shown (see text 
for the definition of references) and agree well with the UltraVISTA measurements within the uncertainties. Right panel; The evolution 
of the number density of galaxies above a fixed mass limit from Ult raVISTA. Open circles denote extrapol ations of the Schechter function 
beyond the data. Black and light blue points are from the NMBS IBrammer etaLl 120111) and PRIMUS HMoustakas et al.|[20]l) surveys 
and agree well with the UltraVISTA measurements. 
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Fig. 8. — Left panel: The evolution of the stellar mass density of star forming (blue) and quiescent (red) galaxies as a function of redshift 
with error bars representing total Icr random errors inclusive of co smic var iance and the errors from photometric uncertainties as derived 
using the MC simulations. At low-redshift the measurements from lBell et al.i (i2003i ) (circle and star) and Bal drv et al.i (i2012i ) (square and 
triangle) are shown. The mass density in quiescent galaxies evolves faster than the mass density in star-forming galaxies, particularly at 
high redshift. Although they dominate the high-mass end of the mass function at z < 2.5, quiescent galaxies do not dominate the overall 
mass density of the universe until z < 0.75 due to a much shallower low-mass-end slope. Right panel: Evolution of the number densities 
of star forming and quiescent galaxies at a fixed mass limit as a function of redshift. 
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0.5) of Figure 6 we plot the SMFs at z 0.1 for the to - 
tal popula t ion fr om th e of studies o f Cole et all ()2001[ ). 
iBell et alJ ([20031) , and iBaldrv et all ([201^. Plotted in 
the middle panel of the fowest-redshift bin are the SMFs 
of star formin g and qui e scent galaxies from IBell et ahl 
(|2003. ) and Baldrv et alJ ()2012[ ). Quahtatively, these are 
similar to our measurements, although we note that the 
selection of star-forming and quiescent galaxies is done 
differently than the UVJ selection in UltraVISTA. 

4.2. The Stellar Mass Density and Number Density 

In the left panel of Figure 7 we plot the integrated stel- 
lar mass density of all galaxies as a function of redshift. 
For consistency with other studies in the literature, the 
stellar mass densities have been calculated by integrating 
the maximum-likelihood Schechter function fits down to 
a limit oi Log{Mstar /^q) = 8.0 at each redshift. We per- 
form this integration using the full maximum-likelihood 
fit, even though at z > 2 the a is not well-constrained. 
In order to account for possible systematic errors caused 
by an underestimate of a due to the limited data depth, 
at all redshifts we also include the uncertainties from the 
maximum likelihood fits with a — -1.2. Therefore, the 
quoted uncertainties in Figure 7, and all subsequent fig- 
ures that use integration of the SMFs span the full range 
of uncertainties for both the a-free and a-fixed SMFs. 

Overplotted in Figure 7 are measurements of the 
stellar mas s density a t vario us redsh ifts from p r evious 
studies bv ICole et all (l200lh (COl), " " 
(B03), iDickinson et"al1 12001 fp 03V 



Bell et al.l (12001 



Glazebrook et al 



(2004) (G04), Dro rv et all l2005l) ( D05]. IRudnick et~ 



(ROej, £bntana et al.' ^2(30^)_(F06 
(E08), [Percz-Gonzale z et al.l 



. , .Eisner et al 
20081) (POSy , 



iMarchcsini ct al.' (2009) (M09), Marchesini et al.' (2010) 
(MIO), Baldry ct al. (2012) (B12), and Moustakas ct alj 
(j^l3) (M13). The substantially larger volume covered 
by UltraVISTA allows for an impressive improvement in 
the uncertainties in the evolution of the stellar mass den- 
sities. Within the uncertainties, the previous measure- 
ments do agree well with the UltraVISTA measurements. 

Stellar mass densities and their uncertainties for star- 
forming and quiescent galaxies have also been computed 
using the same integration method as for the total stellar 
mass density. These are plotted in the left panel of Fig- 
ure 8 as a function of redshift. In general, it is clear that 
at z < 3.5, the stellar mass density of quiescent galaxies 
grows faster than that of star-forming galaxies. We ex- 
plore the implications of this further in § 5. All stellar 
mass densities are listed in Table 2. 

In the right panel of Figure 7 we plot the evolution 
of the integrated number densities of galaxies calcu- 
lated using the same Schechter function fits. The num- 
ber densities are determined down to limiting masses 
of Log(MstaJMQ) = 11.5, 11.0, 10.0, and 8.0, and 
these points are labeled in the figure. Plotted as black 
points are the number densities m easured from the 
NMBS survey (|Brammer "eFalllmil) for mass limits of 
Log{Mstar /^q) = 11-0 and 10.0. These points were mea- 
sured using a catalog constructed in a similar way to 
the UltraVISTA catalog and agree well with the number 
densities in UltraVISTA. Also shown in Figure 7 are the 
integrated numbe r densities at z < 1 calc ulated from the 
PRIMUS survey (|Moustakas et al.ll2013D . which covers 
an area ~ 3x larger than UltraVISTA. These are also 
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Fig. 9. — Best-fit Schechter parameters as a function of red- 
shift for combined SMFs determined using the UltraVISTA data 
(red circles). The square s and circle s are the 2=0 mea- 
surements from ICole et al.l ([2001] ). and IBell et a l. (2003). Also 
shown are parameters from previous r est-frame -optical-sclccted 
SMFs fromlPerez-Gonzalez et all 112003 ). Marche sini et al.l ^2009), 
IMarchesini et al.l 1 12010) , and lCaputi et al. (2011). The SMFs show 
little evolution in M*^^^^ with redshift, but an evolution of 2.70j;g42 
dex in <I>* since z = 3.5. It appears that a flattens with increasing 
redshift; however, this is more likely to be a result of depth of the 
data rather than a true flattening (see text). 

consistent with the UltraVISTA measurements. Simi- 
lar integrated number densities for both the star-forming 
and quiescent populations are shown in the right panel 
of Figure 8, and the values are listed in Table 2. 



5. DISCUSSION 

5.1. The Combined Population 

In Figure 9 we plot the evolution of the Schechter pa- 
rameters M*(jj^, $*, and a as a function of redshift for 
the combined population of both star-forming and qui- 
escent galaxies. All points plotted are from the sin- 
gle Schechter function fits except the 0.2 < z < 0.5 
bin where we have plotted the M*^^^ and combined <l'*s 
of the double Schechter function fits. Also plotted are 
the best-fit Schechter param eters from the l ow-re dshift 
SMFs fromColeet^ pOOl , and iBeUeLaT ('2003), and 
the hig h-rcdshif t rest -f rame- o ptical-selected SMFs from 
Perez- Gonzalez ct al. ( 20081 ) IMarchesini et al.l (|2009t ), 
Marchesini et al. (.2010.) , and lCaputi et al.l (|20lH ). 

Comparison of the parameters in our lowest-redshift 
bin (0.2 < z < 0.5) t o those from the ICole et al.l ()200lD 
and IBell et al.l ()2003() parameters shows good agreement 
for both Mstar and <&*. There is some disagreement in 
a, with our data having a steeper low- mass-end than the 
local SMFs. It is unclear why this is, as the UltraVISTA 
data reach a comparable depth in M^tar as the local stud- 



ies. Part of the discre pancy may be because ICole et ahl 
(|2001f) and IBell et al.l ()2003D fit a single Schechter func- 
tion when a double is required. If we compare our double 
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TABLE 1 

Best-fit Schechtbr Function Parameters for the SMFs 
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TABLE 2 

Number and Stellar Mass Densities: 



Redshift Density log (^) > 8 log (^) > 9 log (^) > 10 log (^) > 11 log (^) > 11.5 



0.2 < z < 0.5 


log ('■/a) 


-1.49±0.06 


-1.89±0.07 


-2.35±0.09 


-3.25±0.13 


-^•52l°:i| 




log(»?Q) 


9 iq+0-07 
^■^^-0.06 


-2.45±0.07 


-2.66±0.09 


-3.33±0.13 


-4.601°:?° 




log ('?SP) 


-1.58±0.06 


-2.03±0.07 


-2.65±0.09 


-4 0S+°-13 
^•"°-0.15 


-6^07l°:2i 




log (pa) 


8.41±0.06 


8.40±0.07 


8.35±0.09 


'-■-'^-0.14 


7.08l°:i« 




log (pq) 


8.19±0.06 


8.19±0.07 


8.18±0.09 


7.91±0.13 


7.ooi°:?i 




log (psf) 


7.99±0.06 


7.97±0.07 


7.85±0.09 


7.08t°:l^ 


i_n 91 

O.OU_Q 35 


0.5 < z < 1.0 


log ()?a) 


-1.69±0.04 


-2.00±0.04 


-2.45±0.06 


-3.49±0.11 


-4.88l°:l« 




log(»?Q) 


r, 0^ + 0.26 
-Z.OD_Q 22 


-2.70±0.04 


-2.85±0.06 


-3.60±0.10 


-5.ioi°:l? 




log (vsf) 


'-■'^-0.03 


-2.09±0.04 


-2.66±0.06 


-4 10+°11 
^•-IL'-O.lO 






log (pa) 


8.26±0.03 


8.25±0.04 


8.19±0.06 


7 70+0.11 
'•'•^-0.10 


6.72i°:;9 




log(pQ) 


7.96±0.03 


7.96±0.04 


7.94±0.06 


7.61±0.10 


6.48l°:l? 




log (psf) 


7.97±0.03 


7.95±0.04 


7.84±0.06 


7.06±0.11 


un 94 

0.0»_g j^g 


1.0 < z < 1.5 


log ()?a) 


-2 n4+"-i'^ 


-2 26+°-°^ 


-2.66±0.06 


-3.77±0.12 


-5 4';+0-35 




log (vq) 


.■3 ni +0 05 


-3 ni +0 05 


-3.12±0.06 


-3.93±0.12 


-i^ 74+0.33 
^•'^-0.19 




log(»?SF) 


-1 Q4+010 


9 Qn+0-05 

-^.CSU_Q Q4 


-2.85±0.06 


-4.27±0.12 


fi 3q+0-29 




log (pa) 


8.02±0.03 


8.01±0.04 


7.95±0.06 


7 41+0-13 
'•*^-0.12 


1 Q+0.36 
"■'■^-0.18 




log(pQ) 


7.65±0.03 


7.65±0.04 


7.64±0.06 


7.25±0.12 


c; 03+0.34 

0.00_Q j^g 




log (psf) 


7.78±0.03 


7.76±0.04 


7.66±0.06 




c: 17+0.30 
^■^'-0.19 


1.5 < z < 2.0 


log (va) 


-2 41 +0-37 


9 c;fi+0-19 


-2.88±0.07 


-3.97±0.14 


-5.76l°:3? 




log(»?Q) 


_Q QQ+0.12 
•^••^^-0.04 


-3 40+010 


-3.48±0.07 


-4.29±0.14 


-6.27l°:« 




log(»?SF) 


9 91 +0.21 


-2.52±0.11 


-3.01±0.07 


-4.25±0.14 


-fi 05+°-26 
-D.UO_g 23 




log (pa) 


7 7Q+0-05 
'•'''-0.03 


7 7Q+0.05 
'•"'-0.04 


7.74±0.07 


7 20+°-l^ 
'•^L)-0.13 


c; 01 +0.39 
^■°1-0.21 




log(PQ) 


7.29±0.03 


7.30±0.04 


7.29±0.07 


6.88±0.14 


c 90+0.42 




log (psf) 


7.65±0.04 


7.64±0.05 


7.56±0.07 


6.93±0.14 


c cQ+0.26 
0.00_g 23 


2.0 < z < 2.5 


log (va) 


-3.061°:™ 


oil +0.43 
-•^•-Li_0.11 


-3 30+0-16 


-4.20±0.16 


-5.89l°:3| 




log {nci) 


-3 ^is+o-ss 

-•5-O»_0.41 


-3 fi7+0-22 

•^•^'-o.so 


-3 91+0-11 


-4.82±0.17 


-6^37l°:3J 




log(»?SF) 


Q 1 0+0. 87 

-3.1»_Q Jg 


Q QQ-|-0.52 

-O.ZO_Q j^2 


-3 42+0-18 
•^■^^-0.09 


-4.32±0.16 


-6.041°:^? 




log (pa) 


743+0.11 


7 43+0.11 
'•''•J-O.06 


7 41+0-10 

'■^-'■-0.08 


6.98±0.16 


r: fiq+0.40 
O.Dy_Q 26 




log (pq) 


6.83±0.07 


6.82±0.07 


6.80±0.09 


6.38±0.17 


r: 9-1+0.32 
^■^1--0.41 




log (psf) 


7 01 +0.14 
'■■^-^-0.04 


7 01+0.13 
'••JJ^-0.06 


7 90+0.10 
' ■''='-0.08 


6.86±0.16 


r; r;3+0.46 
0.00_g 27 


2.5 < z < 3.0 


log {va) 


as+0.69 
-^•»»_0.52 


-3^09l°:« 


-3.43t°:l| 


-4.32±0.18 


7+0.28 
-0.0(_(j 27 




log(»?Q) 


-4 02+1-2^ 


-^•05l°:«^ 


-4 19+0-39 


-5.00±0.20 


-6 63+0-45 
-D.DJ_Q 44 




log (r/SF) 


.9 QQ+0-84 
^•'"'-0.53 


-3 20+° '*° 


c;r;+0.21 

-0.00_g j^g 


-4.42±0.18 


-5 fil+0-28 
-O.Dl_g 27 




log (pa) 


7 09+0.13 
' -"'^-0.09 


^ 09+0.13 
' -"'^-o.io 


7.281°:!? 


6.91±0.18 


6.04±0.28 




log (pq) 


6.581°:?? 


6.58l°:f, 


6.571°:?^ 


6.19±0.20 


4 95+° '** 

'i.90_g 4g 




log (psf) 


7 91 +0.16 
'•^-^-0.10 


7 91+0.15 
'•^^-0.11 


7 17+0.13 
'■^'-0.12 


6.82±0.18 


6.00±0.28 


3.0 < z < 4.0 


log (»?a) 




-3 fi^+0-84 


-4 20+°-^® 
^■'^"-0.34 


-5 04+°-22 


-5 80+°-^^ 
-O.OU_g 24 




log(»?Q) 




-5 44+2-19 


-5 46+1-^^ 

a-^D-0.42 


-5 83+°-3^ 


-6 97+0-69 
D-»' -0.49 




log (vsf) 


-3.iill:S 


Q f^7+0.88 
"•^-°'-0.81 


-4 27+°-38 

^■^'-0.35 


^ 11+0.19 
°-ll-0.17 


-5.85±0.24 




log (pa) 




6 fiS+0-31 


f. P.-+0.20 
'-0.15 


6 31+0-21 


5.89±0.26 




log(PQ) 


0.00_Q 22 


5 r;o+0.81 
0.00_Q 22 


r =0+0.62 
0.00_Q 20 


5 44+0-42 
^•^^-0.19 


4 62+°-" 

*-"^-0.53 




log (psf) 


6 59+° ''° 
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Note. — Number densities ;/ are in units of Mpc^3^ while stellar mass densities p are in units of Mq Mpc^"^. Both the 
number and stellar mass densities are calculated for log (Mgtar/MQ) < 13. The listed errors are the total la errors, including 
Poisson uncertainties, the uncertainties from photometric redshift random errors, and cosmic variance. 
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Schechter function fits (ai = -0.53 



+0.16 
0.281 



a2 



-1.37 



+0.01N 

-0.06.' 



to those derived from lBaldrv et al.) ()2012D (ai = -0.35 ± 
0.18, a2 = -1.47 ± 0.05), we find good agreement. 

Figure 9 also shows that there is good agreement be- 
tween our SMFs and previous high-redshift SMFs in the 
literature. There is a significant improvement in the un- 
certainties in the SMFs derived from the UltraVISTA 
catalog, mostly due to the fact it covers an area that 
is a f actor of 8.8 and 11.4 large r th an the area used 
in thelP^rez-Gonzalez et all (|2008[ ) and IMarchesini et al.l 
(j2009D studies, respectively. Figure 9 confirms the re- 
sults of those previous works and shows that within the 
substantially smaller uncertainties, there is still no sig- 
nificant evolution in M*^^^ out to z 3.0. This lack of 
evolution implies that whatever process causes the expo- 
nential tail of the Schechter function, does so in a con- 
sistent way over much of cosmic time. 

Although there is no evolution in M*jjj^, there is 
a substantial evolution in $* from z = 3.5 to z = 
0.0. If we compare to the $* at z =0 from 
iCole et al.l (|2001h . we find that it evolves by 2.bd>tll\ 
dex between z ~ 3.5 and z ^ 0.0. As Fig- 
ure 9 shows, this evolution is stronger than the val- 
ues of 1.22±0.43, 1.76ti ];s2, l-Q^tpjfi, l-89j: i]; ig de x 
measured previo u sly b y jPerez-Gonzalez et al.' ('2008'), 
Marchesini et al.l (l2009f ). IMar chesini et aL, (2010), and 



Caputi et al.l H2011 ). respectively. 



Interestingly, it appears that there is a statistically sig- 
nificant evolution in a up to z = 2, the redshift where 
the data is deep enough that there is still reasonable con- 
straints on a. A flattenin g of the slope with redsh ift was 
also seen in the SMFs of IMarchesini et all (j2009D which 
probe to slightly lower yigtar- Such a flattening in the 
combined population is a natural consequence of the fact 
that there appears to be little evolution in the a of the 
star-forming population, but a flattening in a for the 
quiescent population with increasing redshift (e.g.. Fig- 
ure 5, see § 5.2). UV -selected sam ples suggest a s teep a 
at high-redshift (e.g. IReddy fc St eidcl 2009; Star k et al.l 
I2009t iGonzalez et al.ll2011l:ILee et al..,2012.) . which taken 
at face value disagree with the flattening of the slope 
observed for the K^-selected SMF. As we show in § 5.3, 
UV-selection misses a fraction of the massive galaxy pop- 
ulation due to their quiescence and/or dustiness, so may 
overestima te a for the combine d population. 

Recentlv lSantini et al.l (|2012t ) have measured the faint- 
end slope using ultra-deep HAWK-I Kg data which is a 
better comparison to the UhraVISTA SMFs than UV- 
selected SMFs. They also flnd a steep faint-end slope (a 
= —1.84 ± 0.06) a.t z ^ 2, which taken at face value does 
not agree well with our m easurement. Of cour se, a and 
M*j^^ are correlated, and iSantini et al.l (|2012| ) cover an 
area that is two orders-of-magnitude smaller than Ultra- 
VISTA and therefore have poor constraints on M*^^^^ . At 
z ~ 2, they flnd a value of Log(M*i^^/MQ) = 11.82 ± 
0.28 which almost an order-of- magnitude larger than the 
UltraVISTA value of Log(M*t„JM0) = 10.81±Ej ;;^. This 
makes it is clear th at in order to simultaneously fit the 
ISantini et al.l ()2012D measurements and the UltraVISTA 
measurements, a double Schechter function is required. 

This implies that the decline in a with increasing red- 
shift in the UhraVISTA SMFs is more a refiection of 
the chugging limiting Mstar with increasing redshift. At 



higher-redshift the UltraVISTA Schechter function fits 
become increasingly dominated by the high-mass end 
and does not contrain the double Schechter function 
upturn. Probing further down the low-mass end will 
most likely result in steeper a than have been measured 
with the current depth. We note that given the limit- 
ing Mstar of the UhraVISTA SMFs, and that a is not 
well-constrained at high-redshift we have also performed 
all the fits with fixed a = —1.4, which is much closer 
to the ISantini et all (|2012f ) and UV-selected SMF mea- 
surements. These fixed-a fits are incorporated into the 
uncertainties in quantities such as the stellar mass den- 
sities and so these should still be representative. 

If the a for star forming and quiescent galaxy SMFs 
shown in Figure 5 are representative of the a at lower 
masses than are probed by the current catalog, then it 
suggests that the single Schechter function fit we have 
used will be a poor representation of faint end at in- 
creasing redshift. Because the a of the star forming 
and quiescent SMFs are so different, and the mix of the 
two populations is a strong function of mass (e.g.. Fig- 
ure 6), a double Schechter function parameterization will 
be necessary to describe the faint end of the combined 
population below the mass limit of the current data. If 
so, then the flattening in a with redshift that we mea- 
sure for the combined population is most likely a con- 
sequence of the fact that the mass limit of the survey is 
near the Mstar where the number densities transition be- 
tween being dominated by the quiescent population, to 
being dominated by the star-forming population. A more 
conclusive understanding of the differences in the slope 
of rest-frame-optical-selected and UV-selected SMFs will 
have to wait until deeper rest-frame-optical data is avail- 
able and can select galaxies to comparable limits in stel- 
lar mass as UV-selected samples. 

5.2. The Quiescent Population 

The SMF of quiescent population shows signiflcant 
evolution since z = 3.5. Although there is little change 

-1-0.03 1 on+0.07 j 
0.04' ^■"-'^-0.07' ^^^"^ 



in M*j^^, the $* increases by 0.57 



2.751q 23 dex since z ^ 1.0, 2.0, and 3.5, respectively. 
With tire improved uncertainties from UltraVISTA, there 
is also evidence for a steady increase in the number 
density of even the most massive quiescent galaxies 
(Log(M^tar/Mo > 11.5) since z = 3.5 (Figure 8). The 
uncertainties from earlier studies (e.g., iBrammer et al.l 
1201 It ) were large enough that they could accommodate 
no growth in the number densities of these galaxies since 
z = 2.2. 

At z < 1 the data is complete to a sufficiently small 
Mstar that the upturn in the number density of quiescent 
galaxies at Log(Mstar/Mo) < 9.5 can be clearly seen. 
This upturn is also seen at z = in the quiesc ent popu- 
lation (e.g.. lBen et"alll2003l: iBaldrv et al.ll2"01^ , and the 
UltraVISTA data now conflrms that it persists to at least 
z = 1. In the appendix we show that the existence of the 
upturn is robust to the deflnition of quiescent galaxies, al- 
though its prominence can be increased or decreased de- 
pending on the strictness of the selection. Interestingly, 
the location of the upturn seems to evolve in mass, from 

Log(M,tar/Mo) - 9.2 at Z = 0.75, to Log(Mstar /Mq) - 

9.5 at z = 0.35. 
It has been suggested that the upturn is the result 
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of a population of star-forming satellite galaxies that 
have been quenched in high-d ensity envir o nmen ts (e.g., 
IPeng et all l20Tol[20ll) . The IPeng et al.l ([20lol ) model 
separates the quenching process into two forms: "envi- 
ronmental quenching" , and "mass quenching" , the for- 
mer of which is caused by high-density environments, 
the latter of which is caused by processes internal to the 
galaxy itself. In their model, the Mstar where the up- 
turn occurs is determined by the relative a's of the mass- 
quenched quiescent population, and the environmentally- 
quenched quiescent population (formerly a recently star- 
forming population). In the case that the a of the 
self-quenched quiescent populations and the star-forming 
populations do not evolve with redshift - which is consis- 
tent with our measurement - an evolution in the M^tar of 
the upturn would imply a change in the fraction of galax- 
ies that self-quench as compared to environmentally- 
quench. An evolution to higher masses with decreas- 
ing redshift would imply an increase in the fraction of 
galaxies that are environmental-quenched with decreas- 
ing redshift, i.e., that environmental-quenching becomes 
increasingly more important at lower redshift. 

It is also interesting to examine the evolution of the 
fraction of quiescent galaxies as a function of both Mstar 
and z. These fractions are plotted in the bottom panels 
of Figure 6. At z < 1 we recover the well-known re- 
sult that quiescent galaxies dominate the high-mass-end 
of the SMF (e.g., fe undv ct al. 2006; Il bert et al.l 120101 : 
iPozzetti et~al] l2010i : iBrammer et al.l 1201 ll ). but there- 
after some interesting trends emerge. 

The fraction of quiescent galaxies with LogiMstar/^Q ) 
> 11.0 continues to decline slowly with increasing red- 
shift up to z — 1.5. After z > 1.5 that decline accelerates 
significantly, and we find that by z = 2.5 the fraction of 
quiescent galaxies has decreased to the point where star- 
forming galaxies dominate the SMF at all stellar masses. 
This result is depicted in Figure 10 where we plot the 
Mstar at which quiescent galaxies dominate the SMF (de- 
noted as Mcross) as a ftinction of redshift. Also plotted in 
Figure 10 a re measurements o f M^^n.^.s from other st udies 
(iBell et al. ''2003': 'B undv et al][200alPo^etti et al.l[2010l : 
iBaldrv et al . 2012). These studies use different defini- 
tions of star-forming and quiescent galaxies, so a direct 
comparison to the UltraVISTA measurements is difficult; 
however, we note that they are largely consistent with the 
UltraVISTA measurements in the redshift range where 
the data overlap. 

As Figure 10 shows, the crossing mass at z = 0.35 is 
Log(Msjar/MQ) = 10.55, and the crossing mass evolves 
as Log{Mcross/MQ) oc (l+z)°-^ up to z = 1.5. There- 
after it evolves much more rapidly Log(Mcrosi5/MQ) oc 
(l-|-z)^'^ up to z = 2.5, after which star-forming galaxies 
dominate the full population. 

Despite the sharp change in the evolution of Mcross at 
z = 1.5, the overall increase in the number density of 
the quiescent population (i.e., $*) is fairly smooth with 
redshift. Therefore, the rapid evolution of Mcross at z > 
1.5 is really a reflection of the fact that Mcross moves 
onto the exponential part of the Schechter function at 
high-redshift, whereas it occurs on the power-law part of 
the Schechter function at low redshift. 

Although their fraction diminishes substantially at 
high-rcdshift, it is noteworthy that we find a non-zero 
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Fig. 10. — The Mstar at which quiescent galaxies dominate over 
star-forming galaxies (Mcross) as a function of redshift. Measure- 
ments from other surveys are shown and agree reasonably well with 
the UltraVISTA measurement. Quiescent galaxies dominate the 
high-mass end of the SMF up to z ~ 1.5. Thereafter star-forming 
galaxies quickly become dominant at all Mstar- 
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Fig. 11. — Left Panel: Comparison of the UltraVISTA star form- 
ing SMF at 2 < z < 2.5 and the BM/BX-selected SMF from 
IHeddv fc Steidell II2009I '|. The SMFs show good agreement at the 
high-mass end but the UltraVISTA SMF suggests a shallower low- 
mass-end slope than the BM/BX SMF. Right Panel: Comparison 
of the UltraVISTA star forming SMF at 2.5 < z < 3.0 and the 
LBG-selected SMF from lHeddv fc S teidell ([20Q9D- These show rea- 
sonable agreement at the high-mass end, although the UltraVISTA 
SMF suggests that the LBG selection may miss ^ 50% of the mas- 
sive galaxy population in this redshift range. 

fraction of quiescent galaxies (^ 10-20%) up to z = 3.5. 
A similarl y small, but non -zero fraction was also found 
by Marche smTet al.l (|2010t) in the NMBS. We wiU discuss 
the SEDs of the quiescent population in more detail in a 
future paper (D. Marchesini, in preparation), but we note 
that for some, the best- fit ages are < 1.0 Gyr. If they 
formed most of their stars in a rapid burst, as suggested 
by the best-fit r-model, it suggests that there may be 
a non-zero population of quiescent galaxies that extends 
to as high as z = 6 - 8. Substantially deeper wide-field 
data (or longer-wavelength selection) will be needed to 
find such galaxies at z > 4, as the Kg-band moves blue- 
ward of the Balmer break and quickly becomes inefficient 
at detecting red galaxies. 

5.3. The Star Forming Population 
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Figure 5 shows that in contrast to the SMF of 
quiescent galaxies, the evolution of the SMF of star- 
forming galaxies from z = 3.5 to z = is fairly modest. 
The M*(^,. and a show no significant evolution up to 
z — 3.5, albeit with large uncertainties in the latter. 
Unlike the quiescent population, the number density 
of Log(Mstar/MQ) > 11.5 galaxies shows almost no 
evolution up to z = 3.5 (Figure 8). The only significant 
evolution is in $*, which evolves by 0.45l5J;^i^, l.Oltg J^^, 

and 2.40;l^o'2i "^^x, since z ~ 1.0, 2.0, 3.5. If we compare 
this to the evolution of the quiescent population at the 
same redshifts, we find that the quiescent population 
has grown faster in <J>* by factors of ^ 1.3, 2.4, and 2.2 
since z ~ 1.0, 2.0, and 3.5, respectively. This shows 
that at all redshifts the majority of the growth in the 
combined SMF is due to the increase in the quiescent 
population. 

The non-evolution in M*^^^ and a and the rather 
slow evolution in $* for the star-forming population is 
remarkable if considered in the context of the evolution 
of the star formation rate per unit Mgtar (specific 
star formation rates, hereafter SSFR) over the same 
redshift range. The SSFR of star-forming galaxies with 
hog{Mstar/MQ) = 10.0(11.0) declines by a factor of - 
20(25) since z ~ 2 (e.g., Muzzin et al. 2013), which 
means the growth rate of these galaxies is evolving 
substantially with redshift. That the process of galaxy 
quenching evolves in such a way to keep the shape and 
normalization of the star-forming SMF roughly constant 
over this redshift range, while there is a significant 
decrease in the SSFRs, implies a carefully orchestrated 
balance between galaxy growth and quenching with 
redshift. 

Given that with the UltraVISTA data we now have 
a reasonably well-determined SMF at z = 3.5 for 
star-forming galaxies, it is interesting to compare this to 
measurements of the SMFs of UV-selected star-forming 
galaxies. In principle, our SMF of rest-frame-optical- 
selected star-forming galaxies should be more complete, 
as both UV-bright and UV-faint star-forming galaxies 
will be selected, where the latter are likely to be highly 
dust-obscured galaxies. With few rest-frame-optical- 
selected SMFs for star-forming galaxies at z > 3 in the 
literature, it is still unclear how large the population of 
UV-faint star-forming galaxies is. 

In the left panel of Figure 11 we plot the SMF 
of Kj-selected star-forming galaxies at 2.0 < z < 
2.5, as well as the S MF o f BM/BX selected galaxies 
by Rcddv & Steidel (2009). We have converted the 
iReddv & Stcidcl (2009) SMF from a Salpeter IMF to a 
Kroupa IMF to match UltraVISTA (N. Reddy, private 
communication). The BM/BX galaxies typically spa n 
the redshift range 1.9 < z < 2.7 (jiReddv fc Steidelll2009D . 
which is a reasonable match to the redshift range of the 
K^-selected galaxies. 

Considering the sizable systematic differences that 
can arise from different ways of measuring the SMF 
(see e.g., the inter-c omparison of SMFs in t he literature 
in the appendix of iMarchesini et all [20090 . the SMFs 
do agree reasonably well, particularly at the highest 
masses. This suggests either that the BM/BX selection 
effectively selects the majority of massive star-forming 
galaxies, regardless of their UV-luminosity, or that 
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Fig. 12. — Comparison of the Ks-selected SMF of star-forming 
galaxies at 3.0 < z < 4.0 from UltraVI STA (blue) and other SMF in 
the literature. The IMarchesini et al.' (2009) and Marchesini et ahl 
(2010) SMFs are also Ks-sclectcd samples and agree well with the 
UltraVISTA SMF. The Stark et al. ( 2009,): .Gonzalez et al.. l,201ll '). 
and Lee et al. (2012) SMFs are UV-selected. These agree reason- 
ably well with UltrsVISTA at Log{Mstar /Mq) = 11.0, but it ap- 
pears the UV-selection may miss the most massive galaxies in this 
redshift range. 

massive dusty star-forming galaxies are less abundant 
than UV-bright ones at z ~ 2.3. Qualitatively speaking, 
the UVJ diagram (Figure 3) suggests that the colors 
of star-forming galaxies above the Msta^-completeness 
limit are primarily red as compared to blue, so the 
high completeness of the BM/BX selection for this 
population is somewhat unexpected. 

The low-mass-end of the SMFs do not agree well, with 
the BM/BX population being more abundant than the 
rest-frame-optical-selected population. This is difficult 
to reconcile given that we have adopted strict Mstar- 
completeness limits and that Ks-selection should be 
more robust than UV-selection. Although we can not be 
conclusive, the discrepancy could arise from systematic 
effects due to the fact that these SMFs are de termined in 
quite different ways. IReddv fc Steidell (|2009D deter mine 
the BM/BX SMF from the UV luminosity function 
(LF) by using a subsample of galaxies for which they 
have spectroscopic redshifts and SED-determined Mstar- 
They convert the UV LF to a SMF assuming that the 
distribution of Mstar in the subsample is representative 
of the entire population. This is different than our 
approach of fitting Zphot and yistar for each g alaxy 
individually. Interestingly, IReddv fc Steidell (j2009l) also 
determine a much steeper faint-end slope (a = -1.73 ± 
0.07) for the SMF than other UV-selected samples which 
measure the SMF with similar Zphot and SED-fitting 
techniques as used for the UltraVISTA catalog (e.g . , 
IStark et al.ll200l IGonzalez et al.llMl ILee et al.il20T2l) . 

In the right panel of Figure 11 we compare the Kg- 
selected star-forming SMF at 2.5 < z < 3.0 to the SMF 
of Lym an Break Galaxies (LBGs) from IReddv fc Steidell 
()2009f) . Down to the mass limit of UltraVISTA, the 
shapes of those SMFs agree reasonably well, although 
the number density of K^-selected star-forming galaxies 
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is a factor of ^ 2 higher than the number density of 
LBGs at Log(Mstar/Mo) = 11.0. Taken at face value, 
it suggests the LBG selection may miss approximately 
half massive gala xies at z > 2.5. A simi lar result was 
also obtained by IMarchesini et al.l f2010), who found 
that only 8/14 galaxies in their mass-complete sample 
(Log(Ms(or/M0) > 11.4) of K^-selected galaxies at 3.0 
< z < 4.0 would be selected with U- and B-dropout 
selection. Us ing t he V VDS spectrosc o pic s ample 
ILe Fevre etHI (I2005D and ICucciati et all ()2012D have 
also found that ~ 50% of the star-forming population 
at z > 3 may obey the LBG selection criteria, although 
we note that that result is based on a comparison of 
luminosity functions, not SMFs. 

In Figure 12 we expand the comparison of the Kg- 
selected SMFs and UV-selected SM F using more recen t 
determinations i n the literat ure from 'Stark et al. (2009); 
IGonzalez efall (|2011D and ILee et al. (2012). We also 
compare to the K^-selected total SMF at 3.0 < z < 4.0 
determine d from the MU S YC /F IREWORKS /FIRES 
surveys (IMarchesini et al.l |2009[) and the NMBS 
(jMarchesini et al.l I2010[ ). In general, most K^-selected 
galaxies at this redshift are star-forming galaxies 
(e.g.. Figure 6), so c ompar ing the total SMFs from 
IMarchesini et al.l (|2009l I2010D to the star-forming SMFs 
is still a reasonable comparison. We note that all of 
the SMFs in Figure 12 have been determined with a 
method similar to the UltraVISTA Kg-selected catalog, 
e.g., Zphot from broadband photometry, and M^tar 
from SED-fitting with similar assumptions about star 
formation histories. These SMFs have a slightly higher 
median redshift than the BM/BX and LBG SMFs, so 
we compare to the highest-redshift SMF in UltraVISTA, 
3.0 < z < 4.0. 

Beginning with the K^-selected samples, we find that 
within the errors, all three Ks-selected SMFs agree 
reasonably well, although the region of comparison is 
limited to fairly high Mgtar- Interestingly, all three 
show little evolution in the number density of the most 
massive galaxies (at Log(Mstar/Mo) > 11.5) from z = 
3.5 to z = 0. The largest discrepancy between the 
Ks-selected SMFs is for UltraVISTA and MUSYC at 
Log(Ms(or/M0) = 11.0, where the number density 
from UltraVISTA is a factor of ~ 3 less abundant in 
such galaxies than MUSYC. Comparison of the 1/Vmax 
points shows that this difference is not significant. 

If confirmed, the result that the abundance of 
Log(Mstar/Mo) > 11.5 star-forming galaxies is un- 
changed since z = 3.5 is quite interesting. We note that 
at present we cannot be 100% confident of how real 
this population is. Examination of their SEDs shows 
that the fits are reasonable; however, they are almost 
all extremely red and dusty. It is also difficult to tell if 
emission from an AGN causes their Mgtar to be inflated. 

At present there are no spectroscopic redshifts for 
such galaxies so we cannot rule out the possibility that 
they are extremel y dusty galaxies from lo wer- z (see also 
the discussion in IMarchesini et al.l I2010D . We explore 
this possibility further in the appendix and find that the 
use of an additional very dusty template in Zphot fitting 
can reduce the number density of this population by 
a factor of ~ 2. Given this, at present it may be best 
to consider their abundance an upper limit to the true 
abundance. 
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Fig. 13. — Evolution of the stellar mass density in the universe 
between z = — 8.5. The SMDs determined from UV-selected 
samples are shown at 2 > 3.5. Below z < 3.5 the Ks-selected 
SMDs from UltraVISTA are shown for the total (black), star form- 
ing ( blue) an d quiescent (red) popu lations. The 2: ^ data 
from ICoIe et al. (2001) (triangle), Belfet^ 1 12003 ) (circles) and 
IBaldrv et'aTT l 2012) (squares) are also shown. The SMD in star- 
forming galaxies from the K^-selected and UV-selected samples 
agrees to within Icr suggesting that UV-selected samples account 
for most of the SMD at 2 > 3.5. The dashed gray curve shows 
a simultaneous fit to the total SMD from UltraVISTA at 2 > 1.5 
and the UV-selected samples, and the dashed maroon shows a fit 
to just the UV-selected samples, both of which agree well. 

Returning to the UV-selected samples, comparison 
of the Ks-selected SMFs to the UV-selected SMFs 
also shows reasonable agreement in the limited mass 
range where the surveys overlap. At Log(M gtar/Mf7)) = 
11.0 t he Kg-selected SMF lies between the iStark et al.l 
(|2009f) and ILee et al.l p012.) SMFs. Extrapolation of 
the best-fit Schechter functions of both of those samples 
does not agree well with the abundance of massive 
star-forming galaxies in Ultra VISTA. A s imult aneous 
Schechter function fit to the ILee et al] (|2012D data 
and the UltraVISTA data does not produce a good 
fit. There seem to be two possible reasons for this. 
Either UV-selection misses most of the most massive 
star-forming galaxies at z > 3 due to dust obscuration, 
or that the abundance of these galaxies is overestimated 
in the K^-selected sample due to a very dusty population 
at lower redshift. 

Overall, the SMFs of rest-frame-optical-selected 
samples and UV-selected samples compare reasonably 
well in the mass range where they overlap. There is 
some evidence that the UV-selection is not complete 
for galaxies with Log{Mstar/M.Q) > 11.0; however, 
spectroscopic verification of the massive population 
Kg-selected population will be important to verify this 
result. 

5.4. Evolution of the Stellar Mass Density 

The left panel of Figure 7 shows that the measured 
evolution of the stellar mass density (hereafter, SMD) 
for the combined population is well-determined and in 
excellent agreement with previous determinations. Our 
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Fig. 14. — Left panel: Average stellar mass of galaxies chosen at a fixed cumulative number density from the Schechter function fits. 
The cumulative number densities are chosen so that they correspond to galaxies with hogCMgtar /^q) = 11.5, 11.0, 10.5, and 10.0 in the 
lowest-redshift SMF. The shaded blue region represent regions that require extrapolating the Schechter function beyond the depth of the 
data. Right panel: Derivatives of the growth in stellar mass as a function of redshift. The derivatives separate into a sequence showing 
that the rate mass growth is always larger for lower-mass galaxies at 2 < 2. 

data show that the SMD of the universe was only 50%, 
10% and 1% of its current value at z ^ 1.0, 2.0, and 3.5, 
respectively. 

The left panel of Figure 8 shows the evolution of the 
SMD of the quiescent and star-forming populations sep- 
arately. The SMD in these populations evolves quite 
differently since z = 3.5, a conclusion that was already 
apparent from the comparison of the SMFs themselves. 
The SMD in star-forming galaxies increases at a rate of 
Pstar oc (1 + z)^^-^^° '^, whereas the SMD in quiescent 
galaxies evolves much faster, as Pstar oc (1 -I- z)~''-^*°-^. 
This strong differential evolution in the SMD can also 
be seen in the inset panel of Figure 8 where we plot the 
fraction of the total SMD in star-forming and quiescent 
galaxies as a function of redshift. 

As Figure 7 shows, at z = 3.5 the universe contained 
only ~ 1% of the total Mgtar formed by z = 0. The in- 
set of Figure 8 shows that at that time the fraction of 
quiescent galaxies was small and approximately 90% of 
the total SMD was contained within star-forming galax- 
ies. Since z = 3.5, the fraction of the SMD in quiescent 
galaxies has grown continuously and at z 0.75 the SMD 
in quiescent galaxies became approximately equal to the 
SMD in star-forming galaxies. Perhaps coincidently, this 
equality in SMD between the two types also occurs pre- 
cisely at the redshift when ^ 50% of the total SMD of the 
universe has formed. Thereafter, the SMD in quiescent 
galaxies exceeds the SMD in star-forming galaxies, al- 
though we note that the details of this statement depend 
on the low-redshift comparison samples, as there are no- 
table differ ences between the l ocal studies of Be ll et al.l 
(|2003) and lBaldrv eTall ([20121 ). 

Due to the superior area and depth of the UltraV- 
ISTA data, we can also compare SMDs determined from 
K^-selected samples at z = 3.5 to those determined 
from UV-selected samples at the same redshift. In Fig- 
ure 13 we plot the SMD from the Ultra VISTA data 
as we ll as from the UV-se le cted samples of iStark et all 
(l2009l). iLabbe et all (|2010l ). iGonzalez et al.l (20111 ). and 
iLee et all (I 2012D. We h a ve als o included t he SMDs at 
z = from iCole et"aI1 (pOOl : iBell et al.l (2003.1 and 



iBaldrv et all (pOll . which means that the measure- 
ments of the SMD in Figure 13 span an impressive red- 
shift baseline of z = 0.0 - 8.5. 

If we compare the Kg-selected SMD at z = 3.5 to 
the SMDs o f the UV-se l ected S MDs at z ~ 3. 7 deter - 
min ed from IStark et al.l ()2009[ ). iGonzalezeiin (|201lD . 
an d ILee et al.l mH ) there is reasonable agreement. The 
SMD from the UV-selected samples agree well with each 
other and are systematically ~ 0.5 dex higher than the 
Kg-selected measurement. The uncertainties in the Kg- 
selected SMD at z = 3.5 are quite large because the data 
just reaches M*^^^ and requires a substantial extrapo- 
lation of the Schechter function. Although lower, the 
Ks-selected SMD does agree with the UV-selected SMD 
within the la uncertainties. 

The reasonable agreement between the total SMD of 
the Ks-selected sample and the UV-selected sample illus- 
trates a key point, namely that at z > 3.5, UV-selected 
samples do select samples that account for most of the 
SMD in the universe. Even though direct comparison 
of the SMFs in § 5.3 showed that UV-selection misses 
approximately half of the massive star-forming galaxies 
and all of the massive quiescent galaxies. Figure 13 shows 
that the SMD in such galaxies at z = 3.5 is fairly small. 
Therefore, UV selection appears to be quite complete for 
the majority of star-forming galaxies at these redshifts. 

Given that our results suggest that UV-selected sam- 
ples should be representative of the SMD at z > 3.5, 
it seems reasonable to fit the SMD over a large red- 
shift baseline. The star formation rate density shows 
a clear decline at z < 1.5 (e.g.. [Hopkins & Bcacom 200^ 
iBouwens et "all 120111 : iSobral et alTl2013t ). so we fit the 
data down to that redshift. Including the UV-selected 
samples with UltraVISTA we find that the total SMD 
evolves as Pstar oc (1 -I- z)~'^-''^'^-'^ from z = 8.5 - 1.5. This 
fit is also in good agreement with a fit to just the UV- 
selected samples which evolve as Pstar oc (1 -I- z)~^ '**°'^. 

The exponent in the overall fit is tantalizingly similar 
to the volume growth of the universe. Why the density of 
stars in galaxies increases at such a rate may be purely 
coincidental, but obviously needs to be investigated in 
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more detail and in the context of models of galaxy evo- 
lution. 

5.5. The Average Mass Growth of Galaxies 

The evolution of the SMFs and SMDs illustrate the 
cosmological evolution of the distribution of galaxies as 
a function of Mstar, and the integrated Mstar within the 
overall galaxy population. While important quantities, 
from the standpoint of modeling the process of galaxy 
evolution, measurements of how individual galaxies as- 
semble their mass may be more useful. Such measure- 
ments are difficult, as they requires the non-trivial task 
of linking galaxies and their descendants through cosmic 
time. 

Several studies have argued that one approach for 
linking galaxies to their descendants is to select 
galaxies at a fixed constant number densi ty (e.g., 
Ivan Dokkum et all 120101: iPapoyich et al.l [201ll ). More 
recentlv, iBrammer et ahl (12011) argued that selecting 
galaxies at a fixed cumulative number density is a bet- 
ter approach, as it is single- va lued at all Mgtar ■ In a 
recent paper, iLeia et all (|2013| ) tested this method on 
semi-analytic models (SAM) of galaxy formation. They 
found that selection at fixed cumulative number density 
recovered the mass evolution of the population to an ac- 
curacy ~ 0.15 dex. They found that the fixed cumu- 
lative number density approach may underpredict the 
mass growth of high-mass galaxies, and overpredict the 
mass growth of lower-mass galaxies, although they note 
that this result depends on the evolution of galaxies in 
the SAM model, which does not properly reproduce the 
SMF as a function of redshift. 

Here we perform fixed cumulative number density se- 
lection using the UltraVISTA SMFs in order to measure 
the average mass evolution of the population. We have 
chosen four fixed cumulative number densities to follow. 
These cumulative number densities are chosen to cor- 
respond to the cumulative number density for galaxies 
with Log(Mstar/M0) = 11.5, 11.0, 10.5, and 10.0, in the 
lowest-redshift SMF bin (0.2 < z< 0.5). 

In Figure 14 we plot the Mstar at these four fixed cumu- 
lative number densities out to 2: = 3. The solid shaded 
region in Figure 14 represents the Msjor-completeness 
limits of the survey. Cumulative number densities that 
require extrapolation of the Schechter function to Mstar 
below the completeness limits are shown as open circles. 

The UltraVISTA data is sufficiently complete, and the 
Mstar evolution sufficiently slow that we measure the 
mass growth of the Log{Mstar /Mq) = 11.5 population 
out to z 3.0. This population demonstrates a re- 
markably slow growth in Mstar-, increasing by only by 
0.3 dex since 2: ~ 3, and 0.2 dex sin ce 2: ~ 2. Using the 
same method ' Brammer et al.l ()2011h measured a growth 
of 0.17 dex for this population since z ~ 2, in excellent 
agreement with our finding. 

With the caveat that increasingly larger extrapolations 
are required. Figure 14 suggests that at z < 2 the growth 
of galaxies with lower Mstar is always faster than galaxies 
with higher Mstar- This can be seen clearer in the right 
panel of Figure 14 where we plot the derivative of the 
Mstar-growth curves. These curves present a remarkable 
sequence at all z < 3, where the derivatives are always 
higher for lower-mass galaxies. 

Figure 14 encapsulates what has already been shown 



by many previous studies, namely that there is a "down- 
sizing" of the galaxy population such that the highest- 
mass galaxies assemble most of their Mstar at high- 
redshift, whereas lower-mass galaxies grow more slowly 
over cosmic time. One of the most interesting implica- 
tions of the curves in Figure 14 is that although low-mass 
galaxies grow at higher rates than higher-mass galaxies 
at z < 2, they do not "overtake" their more massive coun- 
terparts. This is simply because the most massive galax- 
ies have assembled such significant amounts of Mstar at 
very early times. It also implies that at redshifts higher 
than probed by the current data, the mass assembly rate 
of the most massive galaxies must be extremely rapid. 

For example, galaxies with Log{Mstar/MQ) = 11.5 at 
z ~ 0.35 have to assemble as much mass at z < 3 as they 
do at z > 3. The amount of cosmic time between < z < 
3 is ~ 5x more than at 3 < z < 00, which implies that 
even if those galaxies begin forming star shortly after the 
big bang, their mass assembly rates must be substantially 
higher in the past, and must be substantially higher than 
those of low-mass galaxies at any epoch. 

6. SUMMARY 

In this paper we have presented the SMFs of star- 
forming and quiescent galaxies to z = 4 using a Kg- 
selected catalog of the COSMOS/UhraVISTA field. The 
catalog is unique in terms of the large areal coverage 
(1.62 deg2) with a significant depth (Ks,tot = 23.4, 90% 
completeness). The high quality of the data allows for 
arguably the best measurement of the SMFs over a large 
redshift baseline to date. 

The total SMFs agree well with previous measurements 
at 0.2 < z < 3.5, particularly at the high-mass end where 
the wide-field UltraVISTA data provides a substantial 
improvement in the statistical uncertainties. We find no 
significant evolution in M*^^,. out to z = 3.5, although 
the uncertainties in M*j^j, are large at z > 2.5. There is 
also a significant evolution in $* out to z = 3.5. These 
results are also consistent with the results of lllbert et al.l 
(j2013), who also computed the SMFs using the Ultra- 
VISTA data. Most of the evolution in the total SMFs 
is driven by the quiescent population, which grows ap- 
proximately twice as fast as the star forming population 
at all redshifts. Integrating the SMFs we find that the 
SMD of the universe was only 50%, 10%, and 1% of its 
current value at z ^ 1.0, 2.0, and 3.5, respectively. 

Classification of star-forming and quiescent galaxies 
was performed using the rest-frame UVJ diagram. The 
SMFs of these populations evolve quite differently out 
to z =4. The quiescent population evolves much 
faster than the star-forming population, and its growth 
drives most of the growth in the combined SMF. The 
SMD contained in star- forming galaxies grows as (1 -I- 
^^-2.3±o.2 -^^hgreas the SMD in quiescent galaxies grows 
much faster, as (1 + z)^^'^^°-^. The fraction of the total 
SMD contained in quiescent galaxies increases with de- 
creasing redshift, and at z ^ 0.75, the SMD in quiescent 
galaxies becomes equal to that in star-forming galaxies. 
This equivalence in SMD occurs at the redshift where ^ 
50% of the current SMD has formed. 

Starting from low-redshift we find that quiescent galax- 
ies dominate the SMF at high-masses, but that domi- 
nance declines to higher masses with increasing redshift. 
At z > 2.5, star-forming galaxies dominate the SMF at 
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all stellar raasses. 

Comparison of the SMFs of the K^-selected star- 
forming galaxies to the SMFs of UV-selected star- 
forming galaxies at2.5<z<4.0 shows reasonable agree- 
ment. It appears that the UV-selected samples do miss 
approximately half of the population of Log{Mstar /^q) 
> 11.0 star-forming galaxies at z > 2.5 due to the fact 
that they are very dusty; however, we note that this pop- 
ulation of massive dusty star- forming galaxies does need 
to be spectroscopically confirmed. 

Comparison of the SMD for the K^-selected and UV- 
selected samples shows that they agree within the un- 
certainties, which implies that even if the massive dusty 
population exists, it contributes relatively little to the 
total SMD at z > 3.5. This suggests that UV-selected 
samples at z > 3.5 are likely to be representative of the 
majority of the SMD in the universe. Given this con- 
sistency we combined the UltraVISTA SMDs with UV- 
selected SMDs from the literature and fit the evolution 
from 1.5 < z < 8.5. We find that the SMD evolves as 
(1 + z)^-^^'^-^, similar to the volume growth of the uni- 
verse. 

We also perform selection at fixed cumulative num- 
ber density to measure the average growth in the Mstar 
of galaxies. We find that at z < 2, the derivatives of 
the Mstar growth are always larger for lower-mass galax- 



ies, which shows that since that time galaxy growth is 
mass-dependent and primarily bottom-up. In the fol- 
lowing appendix we take a closer look at the effects of 
the assumptions made in SED modeling on the SMFs. 
We also examine the effect of the different definitions of 
star-forming and quiescent galaxies on the SMF. 
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APPENDIX 

THE EFFECT OF SPS MODELS, METALLICITY AND STAR FORMATION HISTORY ON THE SMF 

The determination of stellar masses from SED modeling requires making assumptions about various quantities that 
cannot be well-constrained by broadband photometry such as the SFH, dust att enuation law, and metallicity. De- 
pending on the assumptions made, systematically different estimates of Mstar (e.g.. lMaraston et al.ll2006l iMuzzin et al.l 
[2009a, b; Conroy ct al. 2009) and the SMF (e.g., Jvlarchesini et al. 2009) can result for a given dataset. Previous work 
has shown that typi cally the largest systema tic uncertainties in Mstar and the overall SMF arise from the choice of 
the SPS model fe.g.. iMarchesini et al.,.2009.: .Muzzin et al.ll2009af) . with metallicity and the choice of dust attenuation 
law being lesser, although not unimportant s ource s of systematic error. In one of the first works to do a careful study 
of these systematic errors, IMarchesini et al.l (|2009[ ) found that systematic uncertainties in the SMF from the MUSYC 
survey were approximately as large a contribution to the overall error budget as the random errors. Given that the 
UltraVISTA dataset covers an area > 10 times larger than the MUSYC survey, and has superior photometric data, 
in this appendix we re-examine possible sources of systematic uncertainties on the SMF, with the expectation is that 
they may dominate the overall error budget. 

We explore four different possible sources of systematic error. For each case considered we re-perform the SED 
fitting, generating a completely new catalog of Mst ar for each galaxy. The first source of systematic error we test is to 
change the SPS model from BC03 to the models of 'MarastonI (|2005l ). hereafter M05. The M05 models have a different 
treatment of the thermally-pulsating asymptotic horizontal branch stars (TP-AGB) which, amongst other differences 
with BC03 causes them to produce Mstar that are typically a factor of ~ 0.65 lower than those from the BC03 models 
(;e.g..lWuvts eFaIll2007l: rMuzzin et al ] |2009aHMarchesini et al.|[2009l). W hether these models are a better treatment of 
this phase remains an open issue fe.g.. iKriek et al. l l2010HZibetti et al.l[20T3h . 

We also explore the effect of leaving metallicity as a free parameter. The FAST code allows four different metallicities 
for the BC03 models, Z = 0.004, 0.0 08, 0.02, and 0.05, whic h span a range of sub-solar to super-solar. We also explore 
the effect of using a different SFH. Marasto n et al.l (|2010( ) showed that exponential-declining models were likely an 
overestimate of the Mgtar for star-forming galaxies. In order to allow for an increasing SFH, we explore the effect of 
the "Delayed r-model" option in FAST. The Delayed-r is a SFH of the form SFR oc txe"*/'^, which begins with a 
smoother growth in the SFH before the exponential decline. 

Lastly, we explore t he effect of using an additional very dusty template in EAZY when fitting the Zphot- 
IMarchesi ni et al.l (|2010t ) found that with the inclusion of such a template, approximately half of the massive galaxies 
population at z > 3 could be consistent with a somewhat lower ZpUot in the range 2 < z < 3. Given the prevalence of 
such red SEDs at the high-mass end of the SMF at high-redshift, such a template could be important. 

In Figures 15, 16, and 17 we plot the SMFs derived using these different modeling assumptions for the quiescent, 
star forming, and combined galaxy populations, respectively. The panels are organized by row in increasing redshift 
bins, and below the SMFs in each redshift bin we plot the residuals compared to the default SMF. Within those panels 
the shaded region represents the formal uncertainties in the default SMF. The data points are the SMFs determined 
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QUIESCENT GALAXIES: SYSTEMATIC EFFECTS VARYING SED-MODELING ASSUMPTIONS 
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Fig. 15. — First and third rows: Stellar mass functions for quiescent galaxies in different redshift bins. The stellar mass functions have 
been determined using the different SED modeling assumptions in the legend (see text for details). Second and fourth rows: Difference 
in measured number density at a given stellar mass compared to the default SMF. The shaded region represents the formal uncertainty 
(including Poisson noise, cosmic variance, and modeling errors) in the default mass function. The use of the Maraston models and a free 
metallicity affect the quiescent SMFs at levels larger than the formal uncertainties, but the dusty template and different star formation 
history do not. 



with the 1/Ymax method, and the sohd curves represent the best-fit maximum-likelihood Schechter functions. The 
parameters from all SMFs are listed in Table 3. 

Starting with a comparison of the SMFs for the quiescent galaxies in Figure 15, it is clear that overall, the SMFs 
with the different models are reasonably comparable. There are some clear differences for given assumptions and in 
specific ranges of Mstar and z which stand out in the plots of the residuals. Near the high-mass end, the residuals in 
$* can be large, but it is important to note that because we have compared the residuals as a function of $* and not 
Mstan these can appear extremely large at the exponential tail where a small change in Mstar results in a significant 
change in 

The two parameters that have little effect on the quiescent SMFs are the use of the dusty template, and the different 
SFH history. Quiescent galaxies are typically quite old, and mostly dust-free, so it is unsurprising that the precise 
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STAR-FORMING GALAXIES: SYSTEMATIC EFFECTS VARYING SED-MODELING ASSUMPTIONS 
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Fig. 16. — As Figure 15, but for star-forming galaxies. The stellar mass function of star-forming galaxies is not affected by free metallicity 
or a different star formation history. The SMF is affected by the use of the Maraston models at all redshifts, and at the highest-redshift 
the use of the dusty template reduces the number density of massive star-forming galaxies. 

choice of SFH at early times, and the allowance of extra dust are not important in the SED modeling. 

The allowance of metallicity as a free parameter also has little effect at the high-mass end of the SMF, but does 
reduce the number density of very low-mass quiescent galaxies. This has the effect to reduce the significance of the 
upturn in the quiescent SMF seen at the lowest masses in the low-redshift bins. 

As expected, the most significant systematic effect on the SMF comes from using the M05 models instead of the 
BC03 models. At high-redshift the effect is extremely pronounced, and values of Mgtar are typically a factor of ~ 
2 lower. The effect at higher redshift is expected given that quiescent galaxies are younger there, and the TP-AGB 
phase is most apparent when stellar populations are 0.5 - 2.0 Gyr in age (M05). The effect on the SMF from the SPS 
models diminishes with redshift; however, quite surprisingly the M05 models still produce lower number densities at 
the high-mass end at low-redshift. Given that at old ages the BC03 and M05 models are fairly similar, it suggests that 
the best-fit M05 models fit substantially lower ages for massive low-redshift galaxies than the BC03 models. 

In Figure 16 the SMFs for the star forming population are plotted. Figure 16 shows that similar to the quiescent 
galaxies, the choice of SFH has little effect on their SMFs. Allowing metallicity as a free parameter also does not 
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ALL GALAXIES: SYSTEMATIC EFFECTS VARYING SED-MODELING ASSUMPTIONS 
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Fig. 17. — As Figure 15, but for the combined population. Again, the Maraston models and the dusty template have the largest impact 
on the SMFs. 

change the SMFs within the random uncertainties. As expected, including the dusty template does have a noticeable 
effect at the high-mass end of the star forming SMF at high redshift. In both the 2.5 < z < 3.0, and the 3.0 < z < 
4.0 bins, including the dusty templates reduces the num ber density of Log{M star > 11.0 galaxies by 0.2 - 0.4 
dex. This is fully consistent with .Marchesini et al.l (|2010[ ). who found approximate half of their sample of such galaxies 
was consistent with a lower-redshift solution. If these galaxies truly are very dusty galaxies at lower-redshift, it could 
alleviate some of the tension between the measured SMF and the predicted SMFs from models at the high-mass end. 
It is clear that obtaining spectroscopic redshifts and confirmation of the Mstar for this population is an important task 
for future observations. Unfortunately, as discussed in the text, these galaxies are faint and have extremely red SEDs. 
They are usually only well-detected in bands redward of the H-band. Unless they have very strong emission lines that 
are not completely obscured, obtaining redshifts will be challenging. 

In Figure 17, the SMFs for the combined population are plotted. These reflect most of the same trends as were 
already pointed out in the star forming and quiescent SMFs. In summary, we flnd that for the four sources of sys- 
tematic error we have considered, it appears that two have an impact on the SMFs which is larger than the random 
uncertainties, whereas the other two do not. The choice of SFH appears to have little effect on the SMFs of any 
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TABLE 3 

Best-fit Schechter Function Parameters of the SMFs: Different SED-modeling assumptions 
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Note. — This table is available in its full form in the electronic journal. A subsample of the data is shoivn here to 
represent its form. 

type. Furthermore, allowing metallicity as a free parameter does not affect the SMFs in any notable way other than 
to reduce the size of the upturn of the quiescent SMF at very low Mgtar- This has very little effect on the combined 
SMF, because star-forming galaxies dominate at low Mstar- 

The two areas for concern for the SMFs in terms of systematic uncertainties are the potential dusty massive galaxies 
at high-redshift, and the SPS models. Because star-forming galaxies dominate the SMF at all Mstar at z > 2.5, the 
possibility that some of these galaxies are lower- redshift dusty sources has serious consequences of the high-mass end 
of the SMF at high-redshift. Lastly, it has been well-known for several years now that the choice of SPS model is 
the most significan t systemati c unce rtainty in the determination of M^tar and the SMF. Using the best data available 
at the time, Marc hesini et al.l (120091) sho wed that this uncertainty was at least as large as the random uncertainties 
in the data, and .Marchesini et al.l ( 2010[ ) showed that it was the dominant source of error in the NMBS SMF. With 
a substantially larger survey like UltraVISTA, it is clear that we are now completely in the regime where the SPS 
uncertainties dominate the total error budget in the SMFs. Resolving this issue will be critical in order to better test 
models of galaxy formation at all redshifts. 

THE EFFECT OF DIFFERENT SEPARATION BETWEEN STAR FORMING AND QUIESCENT GALAXIES 

ON THE SMF 

Alth ough there are now several determinations o f the SMFs of star for ming and quiescent galaxies in the literature 
(e.g.. Illbert et al.]l2010l : Brammcr et al. 2011; Dommguez Sanchez et all 2011: Mousta kas et al.. 2013t ). the systematic 
uncertainties in these SMFs due to the defi nition of "quiescence" has not been exam ined in detail. Indeed, given 
that galaxies exhibit a range of SSFRs fe.g.. iNoeske et al.|[2007l: iWhitaker et al.ll2012l Muzzin et al. 2013), and that 
there is a "green valley' in the color-magnitude relation, placing galaxies within a binary definition is implicitly an 
oversimplification of the problem. Here we test how varying the definition of star forming and quiescent galaxies affects 
the SMFs of these types. 
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Fig. 18. — UVJ diagram as Figure 3, but with the varied definitions of quiescent galaxies shown as the dotted line. The dotted lines 
represent extreme definitions of star-forming and quiescent galaxies, but therefore bracket the full range of possible stellar mass functions 
for these populations. 



Best-fit S( 'iieciiteii l-'rxcTiox 



TABLE 

'AHA.ME'l 



1 

OF THE SM1-"S; DiEFEHEXT I AM SELECTIOX 



Redshift 


Sample 


Number 


(Mo) 


(10-4 Mpc-=^) 


a 


(10-4 Mpc-3) 


a2 


0.2<z< 0.5 


Q,UVpVJm 


2489 


11.13±0.03 


Q «fi+0-49 


-0.72±0.03 






0.2 < z < 0.5 


Q,UVmVJp 


6769 


11.18±0.03 


12.38l»-™ 


-1.00±0.02 






0.2 < z < 0.5 


Q,UVpVJm 


2489 


10.89±0.01 


17.30±0.b2 


-0.4 






0.2 < z < 0.5 


Q,UVmVJp 


6769 


10.64±0.01 


48.05±0.03 


-0.4 






0.2 < z < 0.5 


Q,UVpVJm 


2489 


10 94+0 05 


14.56+0;g 


-0.33±0.11 


24+0-4'' 


-1 48+0-21 


0.2 < z < 0.5 


Q,UVmVJp 


6769 


10 92+0-05 


23.55 


-0.48±0.11 


1 9q-H.03 


-1.48±0.11 


0.2 < z < 0.5 


SF,UVpVJm 


16057 


10.94±0.03 




-1.33±0.01 






0.2 < z < 0.5 


SF,UVmVJp 


11777 


10.68±0.04 




-1.37±0.02 






0.2 < z < 0.5 


SF.UVpVJm 


16057 


10.89±0.01 


13.5110-15 


-1.3 






0.2 < z < 0.5 


SF,UVmVJp 


11777 


10.58±0.02 


I3.43l°;^t 


-1.3 







Note. — This table is available in its full form in the electronic journal. A subsample of the data is shown here to 
represent its form. "UVpVJm" corresponds to shifting the UVJ selection box by +0.1 mag in U — V and —0.1 mag 
inV — J; "UVmVJp" corresponds to shifting the UVJ selection box by —0.1 mag in U — V and +0.1 mag inV—J. 
The listed errors are the 1<t Poisson errors. 



In order to provide this test we vary the bounds of the quiescent population -within the UVJ diagram by ± 0.1 mag 
in both U - V and V-J. Figure 18 shows an illustration of how the variation of these bounds appear in the UVJ 
diagram. As Figure 18 shows, an alteration of ± 0.1 mag in color is a rather extreme test; however, we have chosen to 
do so in order to see what the maximum systematic uncertainties will be, and how robust the results from the default 
model arc. 

In Figures 19 and 20 we plot the SMFs generated using the new UVJ definition for the quiescent and star- forming 
galaxies respectively. Again, similar to the figures about the assumptions in SED modeling we have included middle 
panels that show the change in $* relative to the default model. All SMFs derived with the different UVJ selection 
are listed in Table 4. 
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Fig. 19. — As Figure 15 but for quiescent galaxies with two different definitions of star-forming and quiescent galaxies. Varying the 
selection has a much stronger affect on the derived mass functions for high masses than low masses because high-mass galaxies are the 
most-quiescent. 

Examining the SMFs for the quiescent population it is clear that the definition of quiescence is quite important 
in the determination of the SMF. Interestingly, it is much more important for galaxies with Log(MstQr/M0) < 11.0 
than galaxies with Log{Mstar/^Q) > H-O. This is because the most massive galaxies are the most unambiguously 
quiescent. Typically they are the ones with the reddest U - V colors within the quiescent box in the UVJ diagram, and 
therefore variations in the box do not change their classification or number density. Figure 19 shows that the change 
in the number density of galaxies with Log{Mstar /^q) < H-O is a clear function of Mstar, and can reach as high as 
0.2 - 0.4 dex at the lowest masses probed. This change is a refiection of the colors of this population, which are bluer 
than the more massive population, and hence their definition as quiescent is more ambiguous. 

Examination of the SMFs for the star forming population with the different UVJ selection (Figure 20) shows 
opposite trends as the quiescent population. The number densities of the lowest-mass galaxies are mostly unchanged, 
whereas the number densities at the highest masses can change by 0.2 - 0.6 dex. This is again a reflection of the 
fact that the colors of the lowest-mass galaxies are typically very blue, making their classification as star-forming 
galaxies unambiguous; whereas the most massive star-forming galaxies have more intermediate colors making their 
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Fig. 20. — As Figure 15 but for star-forming galaxies with two different definitions of star-forming and quiescent galaxies. The effect 
of different definitions of quiescent on the star-forming population is opposite to quiescent galaxies. Low-mass galaxies are not affected 
because they have the highest specific star formation rates, whereas the most massive star-forming galaxies have much lower specific star 
formation rates. 

classification less clear. This comparison of SMFs makes clear that there are issues of interpretation when dividing the 
spectrum of galaxy SSFRs into a binary classification scheme of star forming and quiescent. It shows that comparison 
of the observed SMFs of star forming and quiescent galaxies to those from models of galaxy formation would require 
a careful matching of definitions. 

One result that appears to be robust to the definition of quiescent is the trend of an increasing fraction of quiescent 
galaxies as a function of Mgtar, and a decrease with z. In Figure 21 we again plot the SMFs of quiescent galaxies 
with the different UVJ selection, but show the quiescence fractions in the middle panels. No matter how quiescence is 
defined, the trend that more massive galaxies are more frequently quiescent than lower mass galaxies at z < 2.0 holds. 
Furthermore, the fact that this trend seems to disappear, and our conclusion that star-forming galaxies dominate the 
SMF at all Mstar at z > 2.5 appears to both hold, no matter the definition of quiescence. It appears that at the z > 
2.5, there are very few galaxies consistent with being quiescent, no matter what the definition. 



COMPARISON TO ILBERT ET AL. STELLAR MASS FUNCTIONS 
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Fig. 21. — As Figures 19 and 20, but in the middle panels the fraction of quiescent galaxies as a function of stellar mass is shown. 
Although the absolute fraction of quiescent galaxies changes with the definition of quiescence, the trend that more massive galaxies are 
more frequently quiescent is robust. The result that star-forming galaxies dominate the SMF at all masses at z > 2.5 is also robust to the 
definition of quiescence. 

Recently Ilbert et al. (2013, hereafter 113) have also computed the SMFs of star-forming and quiescent galaxies up 
to z = 4 based on an independently-generated 30 band PSF-matched catalog of the COSMOS/UltraVISTA field. In 
this appendix we make a comparison between those SMFs and the ones derived in this paper. In Figure 22 we plot 
both our SMFs and the 113 SMFs for the total, star-forming, and quiescent populations in the different redshift bins. 
In Figure 23 a comparison of the SMDs derived from the SMFs of the various populations as a function of redshift is 
shown. 

In general there is reasonable agreement between our SMFs and the 113 SMFs for the combined population at all 
redshifts, particularly at the high-mass end. There is some tension on the low-mass end slopes, and this can also be 
seen in the SMDs where the 113 total SMDs are systematically 0.1 - 0.2 dex higher than our derivation, mostly a 
result of a slightly steeper low-mass end in 113, but also partially due to slightly higher overall number densities at all 
masses in several redshift bins. 

Comparing the SMFs of the star-forming and quiescent galaxies between the surveys shows more mixed agreement 
than the total SMFs. In almost all redshift bins the high-mass end of the SMFs of both types do agree reasonably 
well. The exceptions to this are the lowest-redshift bin, 0.2 < z < 0.5, and an intermediate-redshift bin, 1.5 < z < 
2.0. The disagreement at 0.2 < z < 0.5 is surprising and it is not obvious what its origin is. Both catalogs show 
excellent agreement between the Zphot and the zCOSMOS spectroscopic redshifts, which are most complete for high- 
mass galaxies at low redshift. The difference could result from the definitions of a quiescent galaxy; however, it is 
surprising that it should matter as it is at the lowest-redshifts where the definition of a quiescent galaxy is least 
ambiguous. A better understanding of this discrepancy will require an object-by-object comparison between catalogs. 

The agreement at the low-mass end of the star-forming and quiescent SMFs is not as good as the high-mass end, 
with generally the 113 SMFs having shallower a for the quiescent galaxies and steeper a for the star-forming galaxies. 
This is most likely the result of the different definitions of a quiescent galaxy between the studies. 113 define quiescent 
galaxies using a NUV - M(r) vs. M(r) - M(J) color-color space, which is similar, although not identical to the UVJ 
selection that we use. As we showed in the previous appendix, if we adjust the location of the UVJ box to a more 
conservative cut we would produce quiescent SMFs that are in better agreement with those from 113. Likewise, 113 
could most likely accommodate our SMFs if they were to move the location of their UV - optical color box to a less 
conservative cut. 

Comparison of the SMDs in Figure 23 shows that there is reasonable agreement between the surveys for the quiescent 
population at all redshifts. This is because the a is fairly shallow so the SMD is dominated by massive galaxies. For 
star-forming galaxies 113 derive SMDs that are a factor of 2 higher than ours at high-redshift. Again, this is partially 
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Fig . 22. — Comparison of SMFs of the star-forming, quiescent, and combined populations from this paper to those from lllbert e t ahl 
l|2013l ). The shaded regions represent the best-fit maximum-likelihood Schcchter functions and the associated uncert ainties a nd the solid 
points represent the 1/Ymax SMFs. The colored curves in each panel are the best-fit maximum- likelihood SMFs from lllbert"e t al. (2013). 
The combined SMFs show good agreement in most places. There is some disagreement between the SMFs of the star-forming and quiescent 
populations, which is likely from the different definitions for these used (see text). 



due to the steeper a they derive, but is also partially because their overall number densities are slightly higher at 
all yistar- It may not be a surprise that the largest difference is the low-mass-end slopes, as a is always the most 
difficult part of the SMF to constrain. The low-mass-end slope has frequently been a point of controversy in previous 
measurements of the SMFs as it is the location where the photometry is of the poorest quality, hence the Mgtar and 
Zphot are the poorest-constrained. 

Overall, although there are some specific differences, the comparison between the two sets of SMFs and SMDs derived 
from identical data using different methods shows more consistency than discrepancy up to z = 4. This shows that 
there is a reasonable consensus in the SMFs determined with NIR-selected samples, particularly on the high-mass end 
where the S/N of the photometry is highest. The comparison does illustrate two outstanding issues in the accuracy 
of the SMFs that warrant further investigation. Firstly, it is clear that the definition of star-forming and quiescent 
galaxies needs to be defined in a careful and consistent way. A detailed study of both the UVJ diagram and the 
NUV - optical diagram and the locations of galaxies of a given SSFR in those diagrams would be useful for choosing 
boundaries in both that not only correspond better to each other, but also correspond to the best-possible separation 
of star-forming and quiescent galaxies. Secondly, better measurements of a will be important to resolve differences in 
the SMDs. Given that measuring a requires pushing the data to the lowest S/Ns, it is not surprising t hat discrepancies 
between measurements at the low-mass end are common in this type of work (e.g., the discussion in IReddv fc Steidel 
[2009,) . One obvious step forward in measuring a will be the DR2 UltraVISTA data which should offer a substantial 
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Fig . 23. — Comparison of SMDs of the star-forming, quiescent, and combined populations from this paper to those from Hlbert et ahl 
l|2013l ). In general there is reasonable agreement at most redshifts. Notable exceptions are the SMDs in star-forming at high-redshift which 
is primarily due to a steeper a in Ilbert et al. (201^1). There is also some discrepancy in the SMD in quiescent galaxies at low-redshift, and 
the reason for this is not immediately clear given that the data quality is best for bright galaxies at low redshift and both catalogs agree 
very well with the zCOSMOS spectroscopic redshifts. 

improvement in S/N for the faintest sources. 
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